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I .  INTRODUCTION 


In  this  paper  an  update  is  given  to  the  description  of  a 
regenerative  liquid  propellant  gun  code.^  Most  of  the  basic  governing 
equations  for  the  lumped  parameter  model  are  the  same  as  before,  but  are 
repeated  here  for  convenience.  The  derivations  are  in  the  previous 
report.  The  changes  in  the  equations  are  described,  and  the  reasons  for 
the  changes  are  given.  The  code  has  been  used  to  test  the  effect  of 
various  assumptions.  A  number  of  options  have  been  added  to  the  code, 
and  some  options  that  are  no  longer  considered  useful  have  been  removed. 
A  complete  description  of  the  present  group  of  options  is  given.  In  the 
future,  new  options  will  undoubtably  be  added.  The  code  is  compared  to 
experimental  data  for  a  30mm  gun  fixture  at  the  BRL. 

II.  THE  REGENERATIVE  LIQUID  PROPELLANT  GUN 

Fig.  1  shows  a  diagram  of  one  possible  design  of  a  regenerative 
liquid  propellant  gun  with  an  in-line  piston.  The  monopropellant  is 
pumped  into  the  liquid  reservoir  at  the  beginning  of  the  firing  cycle. 

A  primer  is  ignited  and  injects  hot  gas  into  the  combustion  chamber.  As 
the  chamber  is  pressurized,  the  piston  is  pushed  back.  Because  of  the 
piston  area  differential  between  the  two  regions,  the  piston  will  move 
back  even  when  the  liquid  pressure  is  higher  than  the  combustion  chamber 
pressure.  The  pressure  differential  forces  liquid  propellant  through 
the  holes  in  the  piston  face.  The  propellant  ignites  and  burns  in  the 
combustion  chamber.  The  higher  chamber  pressure  accelerates  the  piston, 
leading  to  more  rapid  injection.  Eventually,  the  pressure  pushes  the 
projectile  down  the  gun  tube. 

The  behavior  in  the  combustion  chamber  is  very  complicated.  A 
primer  forces  hot  gas  into  the  chamber,  leading  to  a  gradual  pressure 
rise.  Liquid  jets  are  forced  out  of  the  vent  holes.  These  jets  may 
break  up  into  droplets  because  of  hydrodynamic  forces  or  due  to  impact 
on  the  wall  of  the  chamber.  The  droplets  formed  may  break  up  further  or 
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coalese.  The  propellant  will  eventually  ignite,  and  may  burn  as 
individual  droplets  or  as  envelope  flames.  A  recirculation  flow  will  be 
set  up  in  the  chamber,  further  complicating  the  behavior  of  the  liquid 
jets.  As  the  projectile  starts  to  move  down  the  gun  tube,  gas  (and 
perhaps  some  unburned  propellant)  will  flow  into  the  tube.  Unlike  most 
solid  propellant  guns,  there  is  generally  a  large  abrupt  area  change 
between  the  combustion  chamber  and  the  gun  tube.  This  will  further 
complicate  the  flow  patterns.  These  processes  are  not  well  understood 
even  at  low  pressures. 

The  design  of  the  piston  considered  here  was  used  in  early  test 
2  3 

fixtures,  ’  but  is  not  convenient  for  actual  guns.  The  vent  holes  must 
be  plugged  before  each  shot,  making  rapid  fire  impractical.  Fig.  2 
shows  a  more  practical  gun  design  (concept  VI) .  The  piston  is  wrapped 
around  a  fixed  central  bolt.  The  liquid  propellant  is  injected  through 
the  annulus  between  the  piston  and  the  bolt.  The  present  code  also 
handles  this  type  of  gun. 

Another  design  being  considered  at  the  present  time  allows  the 
central  bolt  to  move  as  well  as  the  piston  (concept  VIC) .  The  code 
cannot  presently  handle  this  case . 

An  alternate  to  the  in-line  piston  configurations  is  the  reverse 
annular  piston  (RAP)  gun.  For  this  case  the  piston  is  wrapped  around 
the  gun  tube  and  moves  in  the  same  direction  as  the  projectile.  Liquid 
is  injected  directly  into  the  gun  tube.  A  separate  code  has  been 
developed  for  this  type  of  fixture,  and  this  will  be  described  in  a 
separate  report. 


III.  BASIC  ASSUMPTIONS 

The  propellant  in  the  liquid  reservoir  is  assumed  to  be  a 
homogeneous,  isothermal  fluid.  Because  of  the  high  pressures  in  a  gun, 
the  liquid  is  considered  to  be  compressible. 
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Figure  1.  A  Regenerative  Liquid  Propellant  Gun  with  Injection  Orifices. 


INNER 


CONTROL 


BEFORE  FIRING 


AFTER  FIRING 

Figure  2.  A  Regenerative  Liquid  Propellant  Gun  with  an  Annular  Piston. 
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The  orifice  flow  through  the  vents  in  the  piston  is  difficult  to 

model.  The  assumption  is  made  that  the  piston  is  infinitely  thin,  so 

the  orifice  behaviour  can  be  ignored.  The  mass  flux  into  the  chamber  is 

normally  computed  assuming  steady  state  Bernoulli  flow.  Comparisons 

with  transient  lumped  parameter,  one -dimensinal ,  and  two-dimensional 

models  of  the  liquid  reservoir  indicate  that  this  is  a  reasonable 
4 

assumption . 

At  present,  there  is  no  reasonable  procedure  for  predicting  the 
flow  patterns  in  the  combustion  chamber.  The  chamber  fluid  is  assumed 
to  be  homogeneous  and  stagnant.  Lacking  any  detailed  information  on 
flow  in  the  chamber,  the  fluid  velocity  in  the  chamber  is  simply  set 
uniformly  equal  to  zero. 

The  propellant  combustion  is  coupled  closely  with  the  fluid  flow  in 
the  chamber,  and  is  also  very  hard  to  model.  As  a  first  approximation, 
the  liquid  is  assumed  to  combust  instantaneously  as  soon  as  it  enters 
the  combustion  chamber,  releasing  all  its  energy. 

There  is  evidence  that  liquid  accumulates  in  the  combustion  chamber 
at  the  beginning  of  the  firing  cycle,  affecting  the  behavior  of  the  gun. 
Later  finite  burning  rate  models  for  the  liquid  propellant  will  be 
introduced,  in  an  attempt  to  model,  at  least  crudely,  the  effects  of 
possible  propellant  accumulation. 

The  primer  is  assumed  to  be  the  same  liquid  as  the  propellant. 
Normally  the  combustion  and  injection  of  the  primer  is  not  considered  in 
detail.  Instead,  the  initial  pressure  in  the  combustion  chamber  is  an 
input  parameter.  Then  the  amount  of  propellant  needed  to  produce  this 
pressure  is  calculated,  and  this  is  taken  as  the  mass  of  the  primer. 

The  piston  and  projectile  are  not  allowed  to  move  during  the  burning  of 
the  primer.  Later  more  complicated  primer  options  will  be  discussed. 
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There  is  usually  a  large  area  change  between  the  combustion  chamber 
and  the  gun  tube .  The  mass  flux  into  the  gun  tube  is  computed  assuming 
steady  state  isentropic  flow.  Since  the  conditions  in  the  chamber  are 
not  well  known,  the  accuracy  of  this  assumption  is  uncertain. 

The  gun  tube  is  also  treated  as  a  lumped  parameter  region,  using 
modifications  of  the  Lagrange  pressure  distribution  that  has 
successfully  modeled  the  behavior  of  solid  propellant  guns'*.  By 
comparisons  with  a  one -dimensional  model,  improved  versions  of  the 
Lagrange  pressure  distribution  for  liquid  propellant  guns  have  been 
developed. ^ 

Below,  we  derive  the  governing  equations  for  the  simple  lumped 
parameter  model  described  above,  assuming  the  gun  of  Fig.  1.  Then  a 
list  of  options  will  be  given,  with  the  appropriate  modifications  to  the 
base  governing  equations . 

IV .  GOVERNING  EQUATIONS 

With  the  above  assumptions,  the  regenerative  gun  behavior  can  be 
modeled  by  13  ordinary  differential  equations.  These  equations  are 
marked  by  numbers  at  the  left  in  square  brackets  ([1]  to  [13]).  Some 
algebraic  equations  are  required  to  compute  the  coefficients  of  the 
ordinary  differential  equations.  The  equations  describe  the  three 
lumped  parameter  regions  (liquid  reservoir,  combustion  chamber,  and  gun 
tube),  the  mass  flux  between  the  regions,  and  the  piston  and  projectile 
motion.  For  historical  reasons,  the  regions  are  numbered  1,  3,  and  4. 
The  RAP  gun  has  an  additional  region,  the  intermediate  chamber,  numbered 
2.  The  derivations  that  have  not  changed  are  in  the  previous  report.^- 
The  notation  is  given  in  the  glossary. 

A.  LIQUID  RESERVOIR 

The  equations  governing  the  piston  motion  are  straightforward.  Let 
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be  the  volume  of  the  propellant  chamber,  and  A^  be  the  area  of  the 

propellant  side  of  the  piston  (this  area  includes  the  vent  holes).  Let 

s__  be  the  piston  travel  (to  the  left)  and  v__  be  the  piston 
t  S  P5 

velocity.  Then 


av^ 

m - -vps  A1  •  a> 

dt 

The  acceleration  of  the  piston  equals  the  force  on  the  piston 
divided  by  the  mass  of  the  piston.  That  is, 

dv  g 

[2]  —  -  —  [p3  (a3  -  Av)  -  (Pl  +  p  )  (Ax  -  Av)]  ,  (2) 

dt  “ps 

where  A3  is  the  area  of  the  combustion  side  of  the  piston,  p3  is  the 
pressure  in  the  liquid  chamber,  p3  is  the  pressure  in  the  combustion 
chamber,  and  Rps  is  the  mass  of  the  piston.  The  quantity  gQ  -  107 
g/s-cm-Mpa  is  a  conversion  constant  to  put  the  acceleration  in  the 

9 

desired  units  of  cm  /s .  By  assumption,  the  pressures  are  constant 

throughout  the  two  chambers.  The  quantity  pps  is  an  empirical 

correction  to  simulate  the  effects  of  frictional  resistance  (see 

options) .  The  piston  is  assumed  initially  to  be  prevented  from  moving 

toward  the  gun  tube.  So  if  the  initial  acceleration  is  negative,  it  is 

changed  to  zero.  Also,  the  quantity  pps  is  solely  a  resistance 

pressure.  If  the  pps  term  is  larger  than  the  other  pressure  terms,  the 

acceleration  is  set  equal  to  zero.  Also,  if  at  a  later  time  the  piston 

reverses  direction,  the  sign  of  p__  is  changed.  The  equation  governing 

ps 

the  piston  travel  is 


ds 


[3] 


ps 


“  V. 


dt 


ps  • 


From  conservations  of  mass, 


(3) 
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(4) 


dV-^ 

dt  dt  V: 

where  p^  is  the  density  of  the  liquid  and  is  the  mass  flux  out  of 
the  reservoir.  The  energy  equation  can  be  rewritten  as 

2 

dpi  Ci  dp^ 

[5]  - ,  (5) 

dt  g0  dt 

where  c^  is  the  speed  of  sound  in  the  liquid. 

Heat  transfer  and  temperature  changes  in  the  liquid  propellant  are 
expected  to  be  small,  and  are  ignored.  Then  the  speed  of  sound  is  given 
by 


Cl  -  J  gQ  K  /  Pl  ,  (6) 

where  K  is  the  adiabatic  bulk  modulus.  For  the  common  liquid 
propellants,  this  can  be  fit  very  accurately  by^ 

K  -  Kx  +  K2  Pl  .  (7) 

The  corresponding  equation  of  state  is 

k2 

o>  “  U  •  <8> 

The  ordinary  differential  equations  [5]  for  the  pressure  is  not  actually 
required.  Once  the  density  is  known,  the  pressure  can  be  computed  from 
the  equation  of  state  (8) . 

Recently  i  '  equation  of  state  has  been  derived  that  also  includes 

g 

the  temperature  dependence.  This  could  be  used  instead  to  include 
temperature  changes  in  the  liquid.  The  new  equation  predicts  that 
temperature  effects  are  minimal  and  predicts  almost  the  same  pressure 


Pl - i(p\/p 


K, 
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dependence  as  eq .  (8). 


The  mass  flux  out  of  the  reservoir  is  assumed  to  be  steady  state 
Bernoulli  flow 

m13  “  CD  J  2S0  Pi  (Pi  *  P3>  •  (9) 

The  discharge  coefficient  is  an  empirical  correction  to  the  equations 
to  take  into  account  frictional  losses.  Since  the  mass  flux  for  a  gun 
firing  is  not  steady  state,  Cp  may  also  be  used  to  approximate  the  time 
delay  in  reaching  steady  state. 

B .  COMBUSTION  CHAMBER 

Now  consider  the  combustion  chamber.  By  analogy  with  eq.  [1], 
dV3 

[6]  - v  A3  ,  (10) 

dt 

where  V3  is  the  volume  of  the  combustion  chamber  and  A3  is  the  area  of 
the  combustion  side  of  the  piston  (including  holes).  Similarly,  from 
conservation  of  mass 

dp3  P3  dV3  mil  -  m3^ 

[7]  - +  -  ,  (11) 

dt  V3  dt  V3 

where  p3  is  the  density  of  the  gas  in  the  combustion  chamber  and  m3^  is 
the  mass  flux  into  the  gun  tube.  The  energy  equation  is 

2 

dP3  c3  dp3  mjj  (hj^3  -  hQ3)  (7  -  1) 

[8]  - +  -  ,  (12) 

dt  gQ  dt  v3  -  b  m3 

where  c3  is  the  speed  of  sound  in  the  chamber,  h^  is  the  enthalpy  of 
the  liquid,  and  hg3  is  the  enthalpy  of  the  combustion  chamber  gas.  This 
is  based  on  the  Noble-Abel  equation  of  state^ 
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P3  “  ^3  ^3  /  (1  ~  b  p3)  ■ 


(13) 


where  T3  is  the  temperature  in  region  3,  Rg  is  the  specific  gas  constant 
(universal  gas  constant  divided  by  molecular  weight  of  the  gas),  and  b 
is  the  covolume.  This  can  also  be  written  as 

P3  -  P3  (7  ~  1)  cv  T3  /  (1  -  b  p3)  ,  (14) 

where  cv  is  the  specific  heat  at  constant  volume  and  7  is  the  specific 
heat  ratio  Cp/cv.  The  corresponding  equation  for  the  speed  of  sound  is 

/e0  ?  P3 

C 3  -  7  -  •  (15) 

P3  (1  -  b  p3) 

Normally,  gas  will  flow  from  the  combustion  chamber  into  the  gun  tube, 
so  the  mass  flux  m-^  does  not  contribute  to  the  energy  equation.  But  if 
the  gas  flows  from  the  tube  into  the  chamber,  there  will  be  an 
additional  source  term.  This  occurs  so  infrequently  and  the  reverse 
mass  flux  is  so  small,  it  is  sufficient  to  just  set  the  mass  flux  equal 
to  zero  in  this  case. 

The  energy  content  of  a  propellant  is  normally  given  in  terms  of 
impetus  A.  Consider  a  quantity  of  propellant  in  a  constant  volume 
closed  chamber.  The  propellant  combusts,  and  the  result  is  a  gas  with 
some  molecular  weight  Mg  and  temperature  TQ  (called  the  isochoric 
temperature) .  The  impetus  is  defined  as 

A  -  Ru  TQ  /  Mg  ,  (16) 

where  R^  is  the  universal  gas  constant.  Impetus  can  be  obtained  from  a 
closed  bomb  experiment,  if  the  proper  corrections  are  made  for  the 
primer  and  the  heat  loss  to  the  chamber  walls.  For  liquid  propellants, 

Q 

the  impetus  is  usually  calculated  by  a  thermodynamics  code.  The  values 
of  the  molecular  weight,  the  specific  heats,  and  the  covolume  of  the 
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product  gases  are  evaluated  at  the  isochoric  temperature  and  then 
assumed  to  be  constant  with  respect  to  temperature  and  pressure.  The 
impetus  does  depend  on  the  initial  loading  density  (grams  of  propellant 
per  unit  volume).  Values  are  chosen  at  a  standard  loading  density  of 
0.2.  The  effects  of  changing  the  loading  density  are  negligable. 

For  our  case,  it  is  more  convenient  to  work  with  the  chemical 
energy 


•l  “  A  /  (7-1)  -  (17) 

The  enthalpy  of  the  liquid  is  given  by 

hLl  “  el  +  Pi  /  ■  (18) 

The  enthalpy  of  the  gas  is  given  by 

hG3  -  cv  T3  +  P3  /  f>3  “  cp  T3  +  b  P3  •  (19) 

C.  GUN  TUBE 

Finally,  the  gun  tube  (region  4)  is  considered.  The  volume  has  the 
standard  equation 

dVA 

[9]  - vpj  ’  (20) 

dt 

where  Vpj  is  the  velocity  of  the  projectile  and  AA  is  the  area  of  the 
gun  tube.  The  rapid  projectile  motion  creates  a  large  pressure 
gradient.  The  standard  approach  is  to  assume  a  Lagrange  pressure 
distribution,^  that  is,  assume  that  the  density  is  constant  with  respect 
to  space.  It  follows  from  the  one -dimensional  continuity  equation  that 
the  gas  velocity  is  a  linear  function  of  distance.  The  gas  velocity  at 
the  base  of  the  projectile  is  equal  to  the  velocity  of  the  projectile. 
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In  solid  propellant  guns,  the  gas  velocity  at  the  other  end  of  the  gun 
(breech)  is  assumed  to  be  zero.  Since  there  is  a  mass  flow  into  the  gun 
tube,  this  is  not  correct  for  our  model.  However,  the  derivation  with  a 
velocity  at  the  entrance  is  much  more  complicated,  and  the  standard 
Lagrange  distribution  is  used  as  a  first  approximation.  Recently  more 
realistic  pressure  distributions  have  been  studied.** 

Integrating  the  momentum  equation, 

M  x2 

p(x)  -  PL  -  (PR  -  Ppj) -  — r  ,  (21) 

2«pj  V 

where  xR  is  the  distance  from  the  tube  entrance  to  the  base  of  the 
projectile,  pL  is  the  pressure  at  the  gun  tube  throat,  pR  is  the 
pressure  at  the  base  of  the  projectile,  M4  is  the  mass  of  the  gas  in  the 
tube,  Mpj  is  the  mass  of  the  projectile,  and  pRj  is  the  resistance 
pressure.  The  latter  is  an  input  parameter,  which  takes  into  account 
the  shot  start  engraving  force  and  the  frictional  forces  between  the 
projectile  and  the  bore.  Evaluating  at  x-xR  yields 

«4 

PR  “  PL - (PR  ’  Ppl)  •  <22> 

2Hpj 

Integrating  from  the  throat  to  the  projectile  base  yields 

*4 

p4  -  p - (pR  .  p  >  ,  (23) 

6MPJ 

where  p4  is  the  space  mean  pressure.  Then  the  above  two  equations  can 
be  solved  for  p^  and  pR. 

Eq.  (22)  and  (23)  are  based  on  the  assumption  that  the  projectile 
has  started  to  move,  creating  a  pressure  differential.  If  the 
projectile  has  not  moved,  then  region  4  is  treated  as  just  an  extension 
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of  region  3,  that  is, 


P4  “  PR  "  PL  “  P3  ' 

Given  the  pressure  on  the  projectile,  the  acceleration  equation  is 

[10]  “  <PR  "  PPJ}  *4  g°  /  **pj  ‘  (25> 

As  with  the  piston,  the  pressure  Ppj  is  only  a  resistive  pressure.  If 
Ppj  PR’  Pr°jectHe  does  not  move.  The  projectile  travel  Spj  is 


given  by 

[11] 

ds  • 

pj 

at  ■  ■ 

(26) 

The  last 

two 'differential  equations  are  exactly  analogous 

to  region  3. 

That  is, 

[12] 

dpA 

ph  dV4  m34 

(27) 

dt 

v4  dt  v4 

where  p 4 

is  the 

space  mean  average  density  of  the  gas  in 

region  4 ,  and 

[13] 

dp4 

c42  dp4  m34  (hG3  -  hG4)  (7  -  1) 

(28) 

dt 

gQ  dt  v4  -  b  m4 

The  average  speed  of  sound  in  the  tube  is 

c4  "  1 

!&o  1  P4 

(29) 

P4  (1  *  b  p4) 

In  early  versions  of  the  code,  the  mass  flux  into  the  gun  tube  was 
approximated  by  steady  state  Bernoulli  flow.  However,  Bernoulli  flow 
assumes  an  incompressible,  isothermal  fluid.  This  is  a  good 
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approximation  for  liquids,  but  not  for  gases.  A  more  reasonable 
assumption  is  isentropic  flow  (that  is,  adiabatic  and  reversible) . ^ 
The  gas  in  the  combustion  chamber  is  considered  to  be  stagnant.  The  gas 
expands  isentropically  into  the  gun  tube  throat.  The  expanded  gas  then 
mixes  with  the  gun  tube  gas  already  present.  The  process  equations  for 
a  Noble-Abel  equation  of  state  are 

T  (1  /p  -  b)^'l)  -  constant  ,  (30) 


and 


p  (1/p  -  b)^  -  constant  .  (31) 

Assume  that  pt  -  pL<  The  pressure  will  equilibrate  much  more  rapidly 
than  the  temperature  or  density.  The  process  equations  (30)  and  (31) 
can  be  used  to  find  the  throat  temperature  Tt  and  density  pt.  Then  the 
one -dimensional  momentum  equation  can  be  integrated  from  the  stagnant 
conditions  in  the  chamber  to  the  throat,  resulting  in 

vt  -  7  2g0  lb(p3-pt)  +  cpTt  [(p3/pt)(7"1)/ir  -  1]  ,  (32) 


and 


m34  “  V  A^  pt  vt  .  (33) 

The  discharge  coefficient  C^'  is  as  before  an  empirical  correction  for 
loss  terms.  Due  to  lack  of  better  information,  this  is  usually  set  to 
one.  If  the  injection  velocity  approaches  the  speed  of  sound,  choking 
may  occur.  Since  all  the  cases  considered  so  far  involve  injection 
velocities  well  below  the  speed  of  sound,  a  test  for  choking  has  not 
been  implemented.  It  is  not  really  known  how  good  an  approximation  eqs. 
(32)  and  (33)  are  to  the  actual  flow  rate  in  a  gun. 
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V .  VENT  OPTIONS 


For  a  lumped  parameter  model,  the  details  of  the  shape  of  the  gun 
system  are  unimportant.  For  an  in-line  gun  system  all  that  is  needed  is 
the  volume  of  the  liquid  reservoir  and  combustion  chamber,  the  liquid 
and  chamber  side  piston  areas,  and  the  vent  area  between  the  regions. 

For  each  option  described  below,  there  is  a  separate  subroutine.  Each 
subroutine  reads  the  input  data  and,  at  each  time  step,  calculates  the 
time  derivatives  of  the  volumes,  the  piston  acceleration,  and  the  vent 
area.  Because  of  this  modular  construction,  additional  options  can  be 
easily  added.  A  number  of  options  have  in  fact  been  added  and  discarded 
since  the  first  report  on  this  code. 

A.  VENT1 


This  option  assumes  the  gun  fixture  of  Fig.  1.  The  equations  have 
been  given  in  the  previous  section.  The  vent  area  is  allowed  to  vary 
over  time.  While  there  is  no  physical  mechanism  that  will  change  the 
vent  area  arbitrarily,  this  capacity  in  the  code  makes  it  easy  to  study 
the  effects  of  making  any  desired  changes  in  vent  area.  The  physically 
reasonable  case  of  constant  vent  area  is  a  special  case  of  this  option. 

A  table  of  vent  areas  is  read  in  as  a  function  of  relative  piston 
travel.  The  maximum  piston  travel  is  calculated  to  match  the  initial 
liquid  reservoir  volume,  and  the  input  table  is  scaled  so  the  last  input 
piston  travel  matches  the  maximum  piston  travel.  The  advantage  is  that 
in  doing  parametric  studies,  the  vent  area  table  can  bo  read  in  using 
either  absolute  or  relative  distances.  Linear  interpolation  is  used  to 
determine  the  vent  area  for  any  given  piston  travel. 

The  vents  may  also  be  considered  as  valves  that  require  a  certain 
pressure  to  open.  If  the  liquid  pressure  is  below  this  input  opening 
pressure,  the  vent  area  is  zero.  It  would  be  more  reasonable  to  open 
the  vents  when  the  pressure  difference  between  the  liquid  reservoir  and 
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the  chamber  reached  a  specified  level,  but  the  above  has  been 

implemented  to  allow  comparisons  with  a  German  code  with  this  vent 
1 2 

opening  algorithm. 

C .  VENT2 

Recent  regenerative  guns  have  been  designed  with  an  annular  piston 
(see  fig.  2).  The  piston  has  a  circular  hole  in  the  center  which 
surrounds  a  central  rod  or  bolt.  The  propellant  is  stored  between  the 
bolt  and  the  piston.  The  bolt  remains  fixed  as  the  piston  moves.  The 
motion  of  the  piston  increases  the  size  of  the  annular  vent,  and 
propellant  is  forced  between  the  piston  and  the  bolt.  The  back  taper  in 
the  bolt  decelerates  the  piston  at  the  end  of  the  stroke. 

As  before,  a  table  is  read  in  of  vent  area  versus  piston  travel. 

The  radius  of  the  central  bolt  is  assumed  to  vary  linearly  between  table 
entries.  The  piston  is  assumed  to  be  infinitely  thin  (the  shape  is 
unimportant  for  a  lumped  parameter  model).  The  area  Ah  of  the  hole  in 
the  center  of  the  piston  is  read  in.  These  guns  often  have  a  grease 
dyke  between  the  piston  and  the  wall.  The  piston  only  touches  the 
chamber  walls  at  the  very  front.  The  pressure  in  the  grease  dyke  is 
very  similar  to  the  chamber  pressure.  The  area  A  of  this  grease  dyke 

O 

is  also  read  in  (if  there  is  no  grease  dyke,  A  -0) . 

6 

From  this  data,  the  radii  of  the  central  bolt  and  the  volume  of  the 
reservoir  can  be  computed.  As  before,  the  table  is  scaled  to  match  the 
initial  reservoir  volume.  Parametric  studies  varying  the  bolt  shape  but 
keeping  the  charge  weight  constant  can  then  be  easily  performed. 

At  any  given  time  step,  the  piston  travel  is  between  two  input 
distances  ^  and  x^.  The  radius  of  the  central  bolt  is  found  by 
linear  interpolation 


L5 


r  -  r 


i-1 


<sps  '  xi-l> 

+  -  (ri  -  r^) 

<xi  *  xi-l> 


(34) 


The  vent  area  is  then  given  by 


K  ~  Ah  “  * 


(35) 


The  rate  of  change  of  the  volumes  must  include  the  central  bolt.  The 
derivative  of  the  radius  of  the  bolt  under  the  vent  is  required 


dr  ^  (r£  -  rj.j) 


dt 


ps 


(xi  -  x^) 


(36) 


After  some  algebra 


dVl  * 

- vps  A1  +  '  tvps  <r  +  r  ri-l  +  ri-l  > 

dt  K  3 


(37) 


dr 


+  (sps  *  xi-l}  <2r  +  ri-l>  —  J  • 


dt 


Similarly , 


dV3  * 

- -  vps  A3  ’  -  tvps  <r  +  r  ri-l  +  ri-l > 

dt  H  3  v 


(38) 


dr 

+  <sps  ‘  xi-l>  <2r  +  ri-l>  —  1 
v  dt 


The  piston  acceleration  is  given  by 


(2) 


dvps  «o 


dt 


“p- 


[p3  (A3-Ah-Ag)  -(Pi+Ppg)  (ArAh)] 


(39) 


The  acceleration  on  the  piston  no  longer  depends  on  the  vent  area.  The 
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other  governing  equations  are  unchanged. 

VENT 3 . 


The  above  option  is  useful  for  parametric  studies  involving  a 
concept  VI  gun.  If  an  actual  gun  fixture  is  to  be  modeled,  a  slightly 
different  interpretation  is  more  useful. 

The  same  data  as  for  the  VENT2  option  are  read  in.  The  piston 
trrvels  are  now  interpreted  as  absolute  numbers.  The  piston  is  assumed 
to  take  up  a  definite  volume,  although  the  exact  shape  is  still 
unimportant.  The  initial  liquid  volume  has  previously  been  read  in. 

(The  liquid  volume  can  be  measured,  or  computed  using  the  specifications 

for  the  piston  and  bolt) .  The  code  then  computes  how  far  the  piston 

must  move  to  inject  the  appropriate  volume  of  propellant.  This  is  taken 
as  the  maximum  piston  travel.  The  piston  is  assumed  to  hit  the  back 
wall  at  this  point,  and  injection  ceases.  This  is  not  absolutely 
correct,  since  the  piston  will  normally  hit  the  back  wall  while  there  is 

still  a  small  amount  of  liquid  between  the  piston  and  the  bolt,  but  the 

errors  will  be  small,  and  the  behavior  at  the  end  of  the  injection 
stroke  has  very  little  effect  on  performance.  The  governing  equations 
are  now  exactly  the  same  as  in  the  VENT2  option. 

An  additional  option  has  been  added.  Early  models  of  the  concept 
VI  gun,  such  at  the  one  at  BRL,  have  a  movable  block  at  the  end  of  the 
liquid  reservoir  mounted  on  Belleville  springs.  The  injection  area 
between  the  piston  and  the  bolt  is  originally  sealed  by  an  O-xing.  As 
the  combustion  chamber  pressurizes,  the  liquid  pressure  rises.  But 
since  the  liquid  is  almost  incompressible,  it  is  difficult  to  move  the 
piston  enough  to  open  the  vent.  So  the  block  moves  backward  and  allows 
the  piston  to  clear  the  O-ring  and  begin  the  injection  process.  Later 
gun  fixtures  use  a  metal-to-metal  seal  instead  of  the  O-ring  and  do  not 
require  a  movable  block.  But  in  order  to  model  the  BRL  gun  fixture,  a 
Belleville  spring  option  is  included. 
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A  table  is  read  in  of  the  distance  the  block  moves  versus  rhe  force 
exerted  by  the  Belleville  springs.  For  the  BRL  system,  this  has  been 
measured.  The  last  entry  in  the  table  is  the  maximum  distance  the  block 
can  move  before  the  springs  bottom  out.  The  area  and  mass  of  the  block 
are  also  read  in.  Two  additional  ordinary  differential  equations  to 
describe  the  block  motion  are  required 


[Bl] 


and 


dvbk 

dt 


Bo 

-  (Pl  Abk  ‘  fbk> 

Mbk 


(40) 


dsbk 

[B2  ] - vbk  ,  (41) 

dt 

where  sbk  is  the  block  travel,  vbk  is  the  block  velocity,  Abk  is  the 
block  area,  Mbk  is  the  block  mass,  and  fbk  is  the  force  exerted  on  the 
block  by  the  Belleville  springs.  When  the  springs  bottom  out,  the  block 
velocity  is  set  equal  to  zero.  The  only  other  equation  that  needs  to  be 
changed  is  for  the  liquid  reservoir  volume 


av^ 

• • ■  +  vbk  Abk  ■ 
dt 


VI.  PISTON  RESISTANCE  OPTIONS 


(42) 


As  the  piston  moves,  there  is  frictional  resistance  between  the 
piston  and  the  chamber  walls.  There  is  no  good  estimate  for  the  size  of 
this  frictional  resistance.  Because  of  the  large  pressures  involved  and 
the  relatively  slow  speed  of  the  piston,  this  effect  is  not  expected  to 
be  important.  Nevertheless,  piston  resistance  options  are  included  so 
the  effects  of  friction  can  be  studied.  Only  one  option  is  presently  in 
the  code . 
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A.  PIS1 


The  resistance  is  considered  to  be  solely  a  function  of  piston 
travel.  A  table  of  piston  travel  versus  resistive  pressures  is  read  in 
The  table  will  be  normalized  to  the  maximum  piston  travel.  The  piston 
resistive  pressure  PpS  at  any  given  piston  travel  SpS  is  found  by 
interpolation.  The  actual  resistive  force  is  PpS  A^.  An  equivalent 
pressure  is  read  in  rather  than  the  actual  force  to  make  it  easier  to 
compare  the  braking  force  exerted  by  friction  with  the  braking  force 
exerted  by  the  liquid  pressure. 

VII.  PROJECTILE  RESISTANCE  OPTIONS 

Similarly,  as  the  projectile  moves,  there  is  friction  between  the 
projectile  and  the  gun  tube.  More  important  is  the  shot  start  pressure 
The  projectile  is  restrained  from  moving  until  the  pressure  at  its  base 
reaches  an  appropriate  value.  The  choice  of  shot  start  pressure  can 
have  a  large  effect  on  performance  and  on  the  maximum  pressures  in  the 
gun.  At  present  there  is  only  one  option. 

A.  PR0J1 

As  with  the  piston,  a  table  of  projectile  travel  versus  resistive 
pressure  is  read  in.  The  table  is  not  normalized,  since  the  maximum 
projectile  travel  is  always  known  (input  parameter).  The  resistive 
pressure  Ppj  for  any  given  projectile  travel  Spj  is  found  by 
interpolation.  If  spj  is  larger  than  the  last  table  entry  for  travel, 
the  last  table  entry  for  pressure  is  chosen.  The  actual  resistive 
force  is  Ppj  A^ .  The  first  entry  in  the  table  is  the  shot  start 
pressure.  The  projectile  will  not  move  until  the  chamber  pressure 
reaches  the  shot  start  pressure. 
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VIII.  DISCHARGE  COEFFICIENT  -  INTO  CHAMBER 


The  discharge  coefficients  needed  in  equation  (9)  are  input 
parameters . 

A.  DISl 

A  table  of  piston  travel  versus  discharge  coefficient  is  read  in. 
The  table  is  normalized  to  the  maximum  piston  travel.  The  discharge 
coefficient  Cp  for  any  given  piston  travel  SpS  is  found  by 
interpolation. 


IX.  DISCHARGE  COEFFICIENT  -  INTO  GUN  TUBE 

The  discharge  coefficients  needed  in  equation  (33)  are  input 
parameters . 

A.  DISl 

A  table  of  projectile  travel  versus  discharge  coefficient  is  read 
in.  The  discharge  coefficient  Cp'  for  any  given  piston  travel  Spj  is 
found  by  interpolation. 

X.  MASS  FLUX  OPTIONS  -  INTO  CHAMBER 


A.  FLUX1 

Steady  state  Bernoulli  flow  is  assumed.  The  mass  flux  is 
computed  using  equations  (9) .  This  routine  also  reads  in  the  number  of 
vent  holes  (there  is  more  than  one  hole  only  for  the  first  gun  design) . 
The  routine  computes  the  Weber  number,  the  Reynold’s  number,  and  the 
third  Lagrange  number.  These  quantities  are  expected  to  effect  the 
behaviour  of  the  injection.  While  this  code  does  not  use  these  numbers, 
they  are  printed  out  as  additional  information  about  the  injection 
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process . 


B .  FLUX2 

Since  the  gun  environment  is  highly  transient,  steady  state  flow 
may  not  be  a  good  approximation.  A  transient  model  for  injection  has 
been  developed  (gun  design  2).  The  basic  assumption  is  that  the 

space  derivative  of  the  mass  flux  pvA  of  the  fluid  in  the  reservoir  is 
zero.  Then  the  one -dimensional  momentum  equations  can  be  integrated 
from  the  back  wall  to  the  vent  with  the  result 

d(pVA)  22  2  f 

[FI] - [0.5  (v  2 -v3  )  +  g0CD2(p1-p3)/Jdx/A  (43) 

dt 

The  integral  of  the  inverse  of  the  cross  sectional  area  is  approximated 
assuming  a  simplified  piston  shape.  This  additional  ordinary 
differential  equation  is  integrated  to  obtain  the  mass  flux  rate  into 
the  chamber.  This  equation  shows  a  very  rapid  rise  to  steady  state. 

C.  FLUX3 

12 

Heiser  has  used  a  transient  injection  model  more  appropriate  for 
the  first  gun  design.  Only  the  flow  in  the  holes  in  the  piston  is 
considered;  the  liquid  in  the  actual  reservoir  is  considered  a  stagnant 
source.  The  result  is  an  ordinary  differential  equation 

2 

dv3  v3 

[FI]  px  L  -  g0  (px  -  p3)  -  0.5  Pl  — -  .  (44) 

dt  cDz 

The  piston  thickness  L  must  be  read  in.  This  option  was  implemented  to 
allow  comparisons  with  the  code  developed  by  Heiser. 
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XI.  MASS  FLUX  OPTIONS  -  INTO  GUN  TUBE 

A.  FLUX1 . 

Steady  state  isentropic  flow  is  assumed.  The  mass  flux  into  the 
gun  tube  is  computed  ^sing  eq.  (32)  and  (33). 

B.  FLUX2 

Steady  state  isentropic  flow  is  still  assumed.  However,  the 
entrance  velocity  vt  is  considered  known.  This  is  calculated  based  on 
the  assumed  pressure  distribution  in  the  gun  tube  (see  gun  tube  options 
TUBE3  and  TUBE4).  The  governing  equation  (32)  is  solved  for  p^  using 
Newton-Raphson  iteration. 

XII.  GUN  TUBE  PRESSURE  DISTRIBUTION  OPTIONS 

Recently  a  number  of  different  assumptions  for  the  gun  tube 
pressure  distribution  were  implemented  and  tested  against  a  one- 
dimensional  code.  The  results  will  be  published  in  an  upcoming  BRL 
report,^  so  the  details  are  not  included  here.  A  brief  description  of 
the  options  presently  in  the  code  is  included. 

A.  TUBE1 

The  standard  Lagrange  pressure  distribution  is  assumed,  that  is, 
the  density  is  constant  in  the  tube,  which  implies  that  the  velocity  is 
linear.  For  the  purpose  of  computing  the  pressure  distribution,  the 
entrance  velocity  vt  is  set  equal  to  zero.  The  momentum  equation  is 
integrated  to  obtain  a  formula  for  the  pressure  p(x)  in  the  tube.  The 
pressure  p^  at  the  base  of  the  projectile  and  the  pressure  p^  at  the  gun 
tube  throat  are  computed  using  eq.  (22)  and  (23).  Since  vt  will  become 
large,  this  is  not  a  very  good  approximation.  However,  for  low 
performance  guns  with  a  relatively  short  gun  barrel  it  is  still 
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accurate . 

B .  TUBE2 

The  velocity  vt  is  allowed  to  be  non- zero,  but  the  time  derivative 
of  vt  is  set  to  zero.  The  two  equations  for  pressure  are  somewhat  more 
complicated,  and  must  be  solved  simultaneously  with  the  equation  for 
isentropic  flow  (32) .  The  three  equations  are  solved  iteratively  to 
obtain  values  for  vt,  pL,  and  p^. 

Since  vt  does  vary  rapidly  during  the  injection  cycle,  this 
approximation  is  still  not  very  accurate.  However,  it  is  much  more 
reasonable  than  simply  ignoring  the  entrance  velocity,  and  will  predict 
performance  much  more  accurately  for  high  performance  cases. 

C.  TUBE3 

The  term  for  the  time  derivative  of  vt  is  left  in  the  pressure 
equations.  An  additional  ordinary  differential  equation  for  vt  is 
derived  from  the  two  pressure  equations.  The  equation  for  steady  state 
isentropic  flow  is  solved  for  p^  (see  FLUX2 ,  gun  tube  injection 
options) .  Then  p^  and  the  time  derivative  of  vt  can  be  found  from  the 
pressure  equations. 

This  is  quite  accurate  until  the  point  when  the  liquid  propellant 
burns  out  (end  of  stroke).  At  this  point  the  one -dimensional  model 
shows  a  rarefaction  wave  moving  toward  the  projectile.  Both  velocity 
and  pressure  profiles  show  a  sharp  slope  break  at  the  wave  front. 

There  is  an  option  to  model  this  wave.  At  the  end  of  piston 
stroke  the  code  assumes  that  a  rarefaction  wave  begins  to  move  down  the 
gun  tube.  The  velocity  of  the  wave  is  the  speed  of  sound  relative  to 
the  fluid.  The  fluid  velocity  is  assumed  to  be  bi-linear,  with  a  knot 
point  at  the  wave  front.  The  fluid  in  front  of  the  wave  does  not  know 
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that  burnout  has  occured.  So  the  fluid  velocity  at  the  wave  front  is 
found  by  linearly  interpolating  between  the  throat  velocity  at  burnout 
and  the  present  projectile  velocity.  After  much  more  algebra,  the 
desired  quantities  can  still  be  computed.  The  approximation  is  still 
fairly  crude.  But  just  assuming  the  velocity  profile  is  bi- linear  with 
a  knot  point  at  about  the  right  location  gives  very  good  agreement  after 
burnout  with  the  one -dimensional  model. 

D .  TUBE4 

The  primary  error  in  the  above  approximation  is  due  to  the 
assumption  that  the  density  is  constant  in  the  gun  tube.  The  one¬ 
dimensional  model  shows  rather  large  changes  in  density,  especially  as 
the  projectile  moves  further  down  the  gun  tube.  So  the  assumption  is 
made  that  the  density  profile  is  linear  (after  burnout,  bi-linear).  The 
velocity  profile  can  still  consistently  be  taken  to  be  linear.  The 
resulting  model  has  very  complicated  algebra.  However,  the  agreement 
with  the  one-dimensional  model  is  extremely  good. 

XIII.  PRIMER  OPTIONS 

In  present  liquid  propellant  guns,  a  solid  primer  is  used.  The 
primer  injects  hot  gas  into  the  combustion  chamber  to  begin  the  firing 
cycle.  The  exact  details  of  the  primer  model  have  only  a  small  effect 
on  pressure,  so  simplifications  have  been  made.  First,  the  primer  is 
assumed  to  be  the  same  material  as  the  liquid  propellant,  to  save  having 
to  read  In  a  new  set  of  properties  for  the  actual  solid  primer.  The 
three  primer  options  then  offer  varying  levels  of  simplification. 

A.  PRIMl 

It  is  assumed  that  the  primer  completely  combusts  before  the  start 
of  the  calculation,  pressurizing  the  chamber.  The  piston  movement  is 
ignored.  The  code  computes  how  much  propellant  would  be  necessary  to 
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create  the  input  chamber  pressure,  and  records  this  as  the  primer 
weight . 

B .  PRIM2 


The  initial  mass  of  the  primer  is  read  in.  Then  the  primer  is 
distributed  in  the  chamber  as  liquid  droplets.  A  droplet  burning  option 
must  be  chosen  (see  below) .  The  size  of  the  droplets  is  determined  from 
the  drop  diameter  in  the  droplet  burning  option.  This  option  can  be 
used  to  mimic  the  delay  in  time  to  reach  the  primer  pressure. 

C .  PRIM3 


The  primer  is  now  considered  to  be  out  of  the  chamber.  The  mass 
of  the  primer  and  the  injection  time  are  read  in.  The  primer  is  then 
injected  at  a  steady  rate  over  the  desired  injection  time.  The  primer 
is  assumed  to  combust  instantaneously  as  it  enters  the  chamber.  During 
the  injection  process,  the  chamber  density  and  pressure  equations  are 
modified 


dt  v3 


[8] 


dp 3  mn3  (ej_  -  hG3)  (y-D 

-  -  .  .  .  +  _ - 

dt  V3  -  b  M3 


(45) 


(46) 


where  nu^  is  the  mass  flux  of  primer  into  the  chamber.  This  option  was 
added  to  facilitate  comparisons  with  the  code  developed  by  Heiser. 

XIV.  GUN  TUBE  HEAT  LOSS  OPTIONS 


In  gun  systems  there  is  a  large  energy  flux  from  the  hot  gas  to 
the  gun  tube.  This  is  important  for  actual  systems  primarily  because  it 
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causes  erosion  of  the  gun  tube.  The  effect  on  muzzle  velocity  and 
pressure  is  on  the  order  of  a  few  percent.  Heat  loss  is  often  treated 
as  an  adjustable  parameter  for  fine  tuning  a  model  to  agree  with 
experiment . 

A.  HEAT1 

The  heat  loss  to  the  gun  tube  walls  is  ignored.  This  causes  a 
slight  increase  in  muzzle  velocity  and  maximum  pressures. 

B.  HEAT 2 

The  gun  tube  temperature  Tw  and  a  heat  loss  factor  are  read  in. 
The  gun  tube  temperature  is  assumed  to  remain  constant  (infinite  sink). 
On  the  short  time  scale  of  the  ballistic  event,  the  tube  wall  will  heat 
up,  but  treating  the  heat  flux  into  the  solid  makes  the  problem  much 
more  complicated.  Then  the  heat  loss  can  be  represented  as 

Qw  -  4  (T4  -  Tw)  /  d4  ,  (47) 


where  T4  is  the  space  mean  temperature  in  the  tube,  d4  is  the  diameter 
of  the  tube,  is  the  heat  transfer  coefficient,  and  Qw  is  the  heat 
loss.^4  The  variation  of  temperature  in  the  gun  tube  is  ignored.  The 
difficulty  is  in  computing  the  heat  loss  coefficient.  In  the  code  the 
correlation  of  Sieder  and  Tate  for  pipe  flow  is  used.  This  assumes 
fully  developed  flow  in  a  smooth  pipe  with  a  nearly  constant  wall 
temperature.  For  highly  turbulent  flow, 

\  -  (k/d4)  0.026  Re08  Pr1/3  („/mw)°  14  (48) 

where  k  is  the  thermal  conductivity  of  the  gas,  Re  is  the  Reynold's 
number,  Pr  is  the  Prandtl  number,  and  (i  is  the  viscosity.  The 
correlation  is  for  steady  state  flow,  while  the  gun  conditions  are 
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highly  transient.  More  important,  the  correlation  has  been  found 
accurate  for  Reynold's  numbers  up  to  about  10^.  Reynold's  numbers  in 
the  gun  tube  can  be  two  orders  of  magniture  higher. 

If  a  HAN  based  propellant  is  completely  burned,  it  will  form  a 
mixture  of  CC>2,  ^0,  and  Procedures  have  been  developed  to  find  the 

viscosity  and  thermal  conductivity  of  a  mixture  of  gases. ^  The 
procedures  were  developed  for  the  low  density  limit.  Fortunately,  for 
high  temperatures,  the  low  density  limit  is  very  accurate  even  for  very 
high  pressures.^  Fits  were  made  for  the  viscosity  and  thermal 
conductivity  of  the  gases  resulting  from  HAN1845  between  1500K  and 
3500K.  Values  for  other  HAN  based  propellants  are  almost  the  same.  A 
separate  low  temperature  fit  was  made  to  determine  the  viscosity  at  the 
gun  tube  wall  ^w.  The  Reynold's  number  is  given  by 

P  v  d4 

Re - ,  (49) 

4 

where  pv  is  the  space  mean  average  over  the  gun  tube.  The  Prandtl 
number  is  given  by 


Pr 


CP  " 


k 


(50) 


The  equation  for  the  space  mean  pressure  is  modified  by 


dP4  Qw  (7  *  1) 

[13] - ... - 

dt  (1  -  b  p4) 


(51) 


Simulations  run  so  far  indicate  a  heat  loss  of  3%  to  5%  for  120mm  guns, 
and  10%  to  12%  for  30mm  guns.  These  numbers  are  in  the  range  reported 
by  Nordheim,  Soodak,  and  Nordheim^  for  their  analysis  of  experimental 
data,  and  comparable  with  the  heat  loss  normally  used  in  solid 
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propellant  gun  modeling.  In  solid  propellant  gun  modeling,  the  heat 
loss  is  often  modified  to  improve  agreement  between  model  and 
experiment.  The  heat  loss  factor  read  in  above  multiplies  eq.  (51),  and 
can  be  used  to  adjust  the  heat  loss  up  or  down. 

XV.  CHAMBER  HEAT  LOSS  OPTIONS 

It  is  not  known  whether  heat  loss  in  the  combustion  chamber  is 
important.  The  surface  area  of  the  chamber  and  central  bolt  are 
relatively  large,  and  the  gas  velocity  past  the  chamber  walls  is 
unknown.  To  model  the  possible  effects  of  heat  loss  in  the  chamber,  a 
totally  empirical  formula  is  implemented. 

A.  HEAT1 

There  is  no  heat  loss  to  the  combustion  chamber  walls. 

B.  HEAT 2 

The  combustion  chamber  wall  temperature  Tc  is  read  in.  The  heat 
loss  is  assume  to  be  determined  by 

Qc  "  hc  <T3  *  Tc>  /  V3  '  <52> 

Note  that  the  heat  transfer  coefficient  hc  has  different  units  than  h^,. 
Since  th<»  flow  in  the  chamber  is  recirculation  flow  rather  than  pipe 
flow,  a  correlation  based  on  the  diameter  of  the  chamber  is 
inappropriate.  Since  there  is  no  way  to  estimate  hQ  without  knowing  the 
flow  pattern  in  the  chamber,  a  table  of  heat  transfer  coefficient  versus 
piston  travel  is  read  in.  The  heat  transfer  coefficient  for  any  given 
piston  travel  is  found  by  interpolation.  The  equations  for  chamber 
pressure  is  modified  by 
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dp3  Qc  (7  -  1) 

[83  - ... - .  (53) 

dt  (1  -  b  p3) 


The  effect  of  heat  loss  can  now  be  studied,  but  the  size  of  the  heat 
loss  is  totally  arbitrary. 


XVI .  AIR  SHOCK 

The  projectile  will  compress  the  air  in  the  tube  as  it 
accelerates.  This  will  have  a  minor  effect  on  the  velocity  of  the 
projectile . 

A.  SH0CK1 


The  resistance  due  to  the  air  shock  in  front  of  the  projectile  is 
ignored. 


fl.  SH0CK2 


Corner"*  developed  a  formula  for  the  air  shock  based  on  the 
Rankine-Hugoniot  relations.  The  initial  pressure  pa  and  temperature  Tfl 
of  the  air  are  read  in.  The  molecular  weight  of  the  air  and  the  ratio 
of  specific  heats  are  read  in.  The  specific  gas  constant  Ra  is 
calculated  (universal  gas  constant  divided  by  molecular  weight) .  Then 
the  resistance  pressure  pg  is  given  by 


Ps  <7+1)  vpj2 

- 1  +  - + 

Pa  ^  So  Ra  ^a 


(54) 


'Pj 


(So  K  Ta) 


0.5 


(7+1)2  VPj2  0  5 
[7  +  - — — ]°‘5 

16  So  Ra  Ta 


This  new  resistance  pressure  pg  is  just  added  on  to  the  frictional 
resistance  pressure  Ppj . 
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XVII.  BUFFER  OPTIONS 


More  recent  gun  fixtures  use  a  water  buffer  to  slow  up  the  piston 
at  the  end  of  stroke  rather  than  a  taper  on  the  central  bolt.  This 
option  was  implemented,  but  has  not  been  checked  carefully  against  an 
actual  fixture.  Some  changes  may  need  to  be  made  to  the  idealization. 

A.  BUFFI 

There  is  no  buffer. 

B .  BUFF2 


When  the  piston  moves  back,  it  enters  a  chamber  filled  with 
liquid.  The  piston  pressurizes  the  chamber,  which  then  slows  down  the 
piston.  The  buffer  fluid  eventually  is  forced  from  the  chamber  to 
ambient  conditions. 

The  area  of  the  piston  that  enters  the  buffer  is  read  in,  as  well 
as  the  area  of  the  exit  holes  and  the  discharge  coefficient.  The 
initial  volume  and  pressure  of  the  buffer  must  be  entered,  as  well  as 
the  bulk  modulus  coefficients  and  the  density  at  atmosphric  pressure. 
Three  additional  ordinary  differential  equations  are  required  for  the 
buffer  volume,  density,  and  pressure. 

The  buffer  fluid  is  treated  as  an  adiabatic  isothermal  fluid  like 
the  propellant  in  the  reservoir.  The  speed  of  sound  is  determined 
analagously  to  eq .  (6)  and  the  mass  flow  out  of  the  buffer  by  steady 
state  Bernoulli  flow  (eq.  (9)).  The  new  differential  equations  are 

dVbf 

[B1]  - v  A  ,  (55) 

dt  * 
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where  Vb£  is  the  volume  of  the  buffer  and  Abf  is  the  area  of  the  piston 
that  enters  the  buffer  volume.  The  density  equation  is 


dpbf  'bf  dVbf  “bf 

[B2]  - ,  (56) 

dt  vbf  dt  vbf 

where  pb£  is  the  density,  and  m^  is  the  mass  flux  out  of  the  buffer. 
The  pressure  equation  is 


[B3  ] 


dPbf 


d*bf 


dt 


So  dt 


(57) 


where  pb£  is  the  pressure  in  the  buffer  and  cb£  is  the  speed  of  sound. 
The  piston  velocity  equation  is  modified  by 


[2] 


dv. 


ps 


dt 


So 

Pbf  Abf  • 

Ss 


(58) 


XVIII.  REPEAT  OPTIONS 

In  systems  studies,  there  is  often  a  physical  constraint  on  the 
class  of  acceptable  solutions.  Two  common  constraints  are  implemented 
in  the  code.  The  code  will  integrate  repeatedly  until  the  constraints 
are  satisfied. 

A.  REP1 

The  code  integrates  once.  This  is  the  standard  option. 

B.  REP2 

A  solution  is  sought  with  a  given  maximum  liquid  pressure.  The 
liquid  pressure  is  adjusted  by  changing  the  vent  area  (or  areas).  The 


desired  liquid  pressure  and  a  first  guess  for  the  vent  increment  is  read 
in. 


The  problem  is  first  integrated  with  the  input  conditions.  If  the 
maximum  liquid  pressure  is  too  low,  the  first  vent  area  is  incremented 
by  the  input  value .  The  other  vent  areas  in  the  table  are  changed 
proportionally.  If  the  liquid  pressure  is  too  low,  the  vent  areas  are 
decremented.  Once  values  are  obtained  both  above  and  below  the  desired 
maximum  liquid  pressure,  interpolation  is  used  to  obtain  the  new  set  of 
vent  areas.  Convergence  normally  occurs  after  a  small  number  of 
integrations . 

C .  REP 3 

A  solution  is  sought  with  a  given  muzzle  velocity.  The  muzzle 
velocity  is  adjusted  by  changing  the  charge  weight  (initial  reservoir 
volume).  The  desired  muzzle  velocity  and  a  first  guess  for  the 
reservoir  volume  increment  is  read  in.  The  same  iteration  procedure  is 
used  to  find  the  initial  reservoir  volume  that  gives  the  requested 
muzzle  velocity. 


XIX.  CHAMBER  PRESSURE  OPTIONS 

It  is  convenient  to  isolate  part  of  the  model  for  study.  When 
comparing  the  model  against  experimental  data,  small  differences  will 
become  large  quickly.  For  instance,  suppose  that  the  chamber  pressure 
is  slightly  high.  This  will  accelerate  the  piston  more  and  increase  the 
rate  of  injection  of  the  propellant.  The  propellant  burns,  and 
increases  the  pressure  rise  in  the  chamber.  So  to  study  the  liquid 
injection  process,  it  is  better  to  disconnect  the  chamber  pressure  from 
the  injection. 

A.  CHAM1 
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The  chamber  pressure  is  computed  using  the  usual  equations. 


B.  CHAM2 

The  experimental  chamber  pressure  versus  time  is  read  in.  At  each 
time  step,  an  interpolated  experimental  chamber  pressure  overwrites  the 
pressure  computed  using  the  ordinary  differential  equations.  This 
allows  a  study  of  the  injection  process  using  the  experimental  data  as  a 
boundary  condition. 


XX.  DROPLET  BURNING  OPTIONS 

So  far  the  liquid  has  been  assumed  to  combust  instantaneously  upon 
entering  the  combustion  chamber.  There  is  evidence,  however,  that 
liquid  accumulates  in  the  combustion  chamber.  So  simple  rules  for  the 
formation  and  combustion  of  droplets  in  the  combustion  chamber  are 
derived. 

The  behavior  of  the  liquid  as  it  enters  the  combustion  chamber  is 
very  complicated.  A  liquid  jet  is  formed  in  the  chamber.  After  some 
time,  the  jet  will  break  up  into  droplets,  ligaments,  or  slugs.  The  jet 
may  impinge  upon  the  walls  of  the  chamber.  At  present,  there  is  no  way 
of  predicting  the  behaviour  of  the  jet.  In  the  model,  the  liquid  is 
assumed  to  instantaneously  form  droplets  of  fixed  size  as  the  liquid 
enters  the  combustion  chamber.  The  droplets  will  then  combust  at  some 
given  rate.  The  initial  size  of  the  droplets  is  an  input  parameter. 

This  simplified  model  has  no  basis  in  reality,  except  to  allow  the  code 
to  determine  the  effects  of  liquid  accumulation. 

In  the  earlier  version  of  the  code,  there  was  an  option  to  allow 
the  droplets  to  decrease  in  size  as  combustion  proceeded.  But  since  the 
initial  droplet  size  is  arbitrary,  this  did  not  really  improve  the 
model.  To  simplify  the  code,  this  option  was  removed.  The  droplet  size 
is  assumed  to  be  uniform  with  respect  to  space.  However,  the  droplet 
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size  may  vary  over  time. 


The  derivation  of  the  relevant  equations  is  very  similar  to  the 
previous  report.  Let  be  the  liquid  mass  in  region  3,  be  the 
liquid  volume,  be  the  gas  mass,  and  VG3  be  the  gas  volume .  Then 


^L3 


(59) 


and 


MG3 


(60) 


where  p ^  is  the  density  of  the  liquid  and  pQ ^  is  the  density  of  the 
gas.  Define  the  porosity 


Then 


(61) 


M-j 


p3 - (1  *€3>  ^L3  +  £3  pG3 


(62) 


Following  eq.  (6),  the  speed  of  sound  in  the  liquid  is 


CL3  “  •>  ®o  K  /  PL3  ’ 

and  following  eq.  (15)  the  speed  of  sound  in  the  gas  is 


'G3 


J  So 


1  P3 


(64) 


PG3  ^  b  PG3^ 

After  some  manipulation,  the  speed  of  sound  in  the  mixture  is  given 
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by 


-f- 


P 3  e3^pG3  cG32)  +  ^  ‘  f3^/^pL3  CL32^ 


(65) 


The  enthalpy  of  the  liquid  is 


”  ei  +  P3  /  ^L3  1  (66) 

where  e^  is  the  chemical  energy  of  the  liquid,  and  the  enthalpy  of  the 
gas  is 


hG3  ”  cp  T3  +  b  P3  •  <67> 

Let  m3  be  the  rate  at  which  the  liquid  in  region  3  is  combusting  and 
forming  the  final  gas  products. 


The  thirteen  equations  derived  in  Section  III  are  the  same  except 
for  the  energy  equations  [8]  and  [13].  The  equation  from  the  previous 
report  is  written  in  a  more  general  fashion  to  make  it  easier  to  apply 
to  different  gun  systems. 


[8] 


,  2 

dpj  P 3  c3 
dt  So  v3 


dV3 

mL3 

“G3 

dt 

pL3 

^G3 

(68) 


S3  So  (7  *  1) 

2  ^  ’ 
qG3  CG3  (!  *  b  Pg3^ 

where  m^  is  the  rate  of  production  of  liquid:  in  this  case. 


mL3  "  “13  '  m3 


M 


L3 


and  m^3  is  the  rate  of  production  of  gas, 


(69) 
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“G3 


nii  -  m 


*G3 

34  +  “p3 

Mo 


(70) 


and  S3  is  the  rate  of  production  of  mass  times  the  change  in  enthalpy 
S3  —  (hj^j  -  hgj)  +  (h^i  *  h^j)  +  (71) 

“p3  <el  *  hG3^  • 

If  there  is  heat  loss  in  the  chamber,  the  modification  term  from  eq. 
(53)  is  added  on. 

In  the  gun  tube,  the  equations  are  exactly  analogous.  The  only 
difference  is  in  the  coefficients  for  the  pressure  equation 


dP4  Ph  c4  dV4  mL4  “G4 

[13] - { - +  -  +  -  + 

dt  gQ  VA  dt  PlA  pG4 


S4  g0  (7  -  1) 


(72) 


1G4  cG4  (1  -  b  Pq4^ 


In  the  gun  tube 


M 


m 


L4 


m/,  +  m 


L3 


34 


(73) 


M, 


mG4  "  m4  +  m34 


^G3 


(74) 


M, 
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M 


s4  "  m4  (hL4  '  hG4^  +  m34  ^hL3  *  hL4^ 


L3 


(75) 


M, 


m34  (hG3  *  ^G4^  • 

m3 

If  there  is  heat  loss  in  the  gun  tube,  the  modification  term  from  eq. 
(51)  is  added  on. 


To  complete  the  system,  ordinary  differential  equations  are 
required  to  determine  the  liquid  density  and  mass  in  the  regions.  The 
corresponding  gas  quantities  can  be  easily  derived.  The  basic 
assumption  is  that  the  pressure  in  the  liquid  and  the  gas  in  a  given 
region  is  the  same.  Since  the  gas  pressure  is  known,  the  liquid 
equation  of  state  can  be  written  in  differential  form  as 


[14] 

dpL3 

o 

is 

1 

(76) 

and 

dt 

CL3  dt 

[16] 

d^L4 

gG  dP4 

-  - 

(77) 

dt 

CL4  dt 

The  mass 


conservation  equations  for  the  liquid  have  already  been  derived 


[15] 


dM 


L3 


dt 


mL3  1 


and 


(78) 


[17] 


dM 


L4 


-  m 


dt 


L4  • 


(79) 
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To  close  the  system,  the  rate  at  which  the  liquid  droplets  are 
combusting  must  be  known.  The  rate  of  surface  regression  is  assumed  to 
be  of  the  form 

a 

rate  of  surface  regression  -  Ap  .  (80) 

Liquid  propellant  burning  rates  have  been  measured  by  McBratney^ ’ ^ . 

The  rate  of  combustion 

m3  “  PL3  ®  APB  •  (81) 

where  S  is  the  total  surface  area  of  the  droplets  in  a  region.  All  the 
droplets  are  assumed  to  have  a  constant  diameter  d.  So  the  total 
surface  area  in  the  region 

S  -  6  VL3  /d  .  (82) 

and  eq.  (81)  can  be  written  as 

m3  “  ^L3  (®/d)  Ap3B.  (83) 

For  region  4,  the  pressure  is  no  longer  constant  over  the  region, 

but  follows  the  Lagrange  pressure  distribution.  But  this  would  be  quite 
complicated  to  keep  track  of,  and  very  little  liquid  goes  into  region  4 
for  most  problems.  So  the  combustion  rate  is  assumed  to  depend  on  the 
average  pressure  p4,  and 

m4  -  mL4  <6/d>  AP4B-  (84) 

The  droplet  options  can  now  be  described. 

A.  DR0P1 

The  liquid  is  assumed  to  combust  instantaneously  as  soon  as  it 


38 


enters  the  combustion  chamber.  This  is  the  model  described  in  section 
III. 


B.  DR0P2 

The  liquid  enters  the  combustion  chamber  and  instantaneously  froms 
droplets  of  diameter  d.  The  size  of  the  droplets  is  fixed  in  time  and 
space.  As  the  droplets  burn,  the  number  of  droplets  may  decrease,  but 
the  size  of  the  droplets  remains  unchanged. 

The  droplet  diameter  and  the  burning  rate  coefficients  are  read 
in.  There  is  some  evidence  that  the  burning  rate  for  liquid 
monopropellants  has  a  sharp  slope  break  around  100  MPa.  So  the  burning 
rate  is  read  in  as  two  fits,  and  the  pressure  is  read  in  at  which  the 
change  is  made  from  the  low  pressure  rate  to  the  high  pressure  rate. 

C.  DROP 3 

The  droplet  size  is  fixed  in  space,  but  not  in  time.  A  table  of 
droplet  diameters  versus  piston  travel  is  read  in.  As  usual,  the  table 
is  scaled  to  the  maximum  piston  travel. 

While  this  model  is  unrealistic,  it  allows  the  modeling  of  various 
accumulation  rates  at  different  times  in  the  firing  cycle. 

XXI,  MASS  AND  ENERGY  BALANCE 

As  a  check  on  the  equations  derived,  both  mass  and  energy  should 
be  conserved.  The  mass  balance  is  straightforward. 


Mt  -  Mx  +  Ml3  +  Mqj  +  M^  +  Mg4  +  Mp 


(85) 


where  Mp  is  the  mass  of  the 
propellant  and  primer.  Mj  should  be  constant  throughout  the 


primer  and  M^  is  the  total  mass  of  the 
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integration. 


The  energy  balance  is  more  complicated.  The  liquid  is  considered 
to  be  an  isothermal  fluid,  and  its  internal  energy  e-^  is  just  the 
chemical  energy  of  the  propellant.  The  total  energy  in  region  1  is 

.  (86) 

If  the  outside  primer  option  is  chosen  (PRIM3)  then 


EP  “  el  *P  • 

The  internal  energy  of  the  liquid  in  region  3  is  given  by 
EL3  "  el  ML3  • 

and  the  internal  gas  energy  is 


(87) 


(88) 


EG3  "  cv  T3  MG3  1 
Similarly  for  region  4 


(89) 


EL4  “  el  ML4  - 


(90) 


and 


EG4  -  <=v  T4  «G4  •  (91> 

The  kinetic  energy  of  the  piston  is  given  by 

EKps  -  °-5  V  V2  >  ■  <92) 


and  the  kinetic  energy  of  the  projectile  by 


(93) 


“pj  -  0.5  Mpj  vpJ2  /  g0  . 

The  liquid  and  gas  in  region  4  are  moving,  and  are  assumed  to  have 
the  same  velocity.  The  kinetic  energy  of  the  liquid  is 

E*L4  “  °‘5  ML4  vpj2  XR  /(3g0>  .  (94) 

where  x^  is  the  distance  from  the  entrance  to  the  gun  tube  to  the  base 
of  the  projectile.  The  kinetic  energy  of  the  gas  is 

EKg4  -  0.5  ^  vpj2  xR  / (3g0)  •  (95) 

The  fluid  in  the  combustion  chamber  is  considered  stagnant.  When 
the  fluid  enters  the  gun  tube,  it  instantaneously  acquires  kinetic 
energy.  To  keep  the  energy  balanced,  the  enthalpy  of  the  liquid  in  the 
gun  tube  is  redefined  as 

hL4  *  el  +  P4  /  ^L4  +  EKL4  /  M4  >  <96> 

and  the  enthalpy  of  the  gas  in  the  gun  tube  is 

hG4  -  Cp  T4  +  b  p4  +  EKq4  /  M4  .  (97) 

The  heat  loss  to  the  chamber  and  gun  tube  is  energy  that  leaves 
the  system.  The  total  energy  loss  is  approximated  by  finite 
differences.  The  code  is  set  up  to  print  out  information  at  designated 
time  intervals.  After  each  time  interval,  the  heat  loss  from  the 
chamber  is  approximated  by 

EH3(t)  -  EH3(tQ)  +  Qc  V3  (t  -  tc)  ,  (98) 

where  Qc  and  V3  are  averages  of  the  values  at  the  present  time  t  and  the 
old  time  tQ.  Similarly 
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EH4(t)  -  EH4(tc)  +  Qw  V4  (t  -  tQ)  . 


(99) 


The  energy  lost  through  friction  is  approximated  in  a  similar 
fashion.  The  energy  lost  is  equal  to  force  times  distance,  or 

EFps(t)  -  EFps(t0)  +  ppg  A3  (sps(t)  -  sps(tD))  .  (100) 

for  the  piston  and 

EFpj (t)  -  EFpj(t0)  +  Ppj  A4  (spj(t)  -  spJ(t0))  ,  (101) 

for  the  projectile. 

Then  the  total  energy  of  the  system  is 

eT  ”  EL1  +  EP  +  EL3  +  EG3  +  EL4  +  EG4  +  ^s  +  (102) 

EKpj  +  EKLA  +  EKG4  +  EH3  +  EH4  +  EFps  +  EFpj  . 

This  should  be  constant  throughout  the  integration. 

XXII.  NUMERICAL  METHOD 


Because  of  the  length  (over  100  pages)  the  code  listing  is  not 
given.  Appendix  A  gives  comments  from  the  beginning  of  the  code.  This 
is  a  complete  description  of  the  input  options  and  the  notation  actually 
used  in  the  code  for  the  input  variables.  A  copy  of  the  code  is 
available  from  the  author  on  request. 

The  ordinary  differential  equations  derived  above  are  solved  using 
EPISODE.1^  This  is  a  robust  and  efficient  code  for  the  solution  of 
ordinary  differential  equations.  EPISODE  proper  consists  of  the 
subroutines  DRIVE  to  SING.  A  few  minor  changes  have  been  made  in  in  the 
subroutines  DRIVE  and  TSTEP. 
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The  subroutine  DRIVE  controls  the  integrator.  The  error  control 
is  defined  in  DRIVE,  based  on  the  input  error  bound  EPS.  Originally  the 
error  control  could  be  either  relative  or  absolute.  A  relative  error 
control  can  be  wasteful  if  some  of  the  quantities  being  integrated 
become  negligably  small.  However,  an  absolute  error  control  is 
inaccurate  if  the  quantities  integrated  vary  widely  in  magnitude.  A 
solution  to  this  problem  was  developed  during  work  on  integrating 
chemical  kinetics  networks,  where  the  concentrations  of  species  differ 
by  many  orders  of  magnitude.  This  is  a  semi -relative  error  control. 

If  the  quantities  are  above  some  cutoff  value  (SREC)  a  relative  error 
control  is  used.  If  they  are  below  SREC,  an  absolute  error  control  is 
used.  While  an  absolute  error  control  rould  be  adequate  for  the  simpler 
equations  now  being  considered,  the  semi-relative  error  control  has  been 
left  in  the  code.  SREC  has  normally  been  chosen  three  orders  of 
magnitude  smaller  than  EPS. 

Also  in  DRIVE  are  various  error  messages.  Among  other  things,  the 
code  prints  out  a  warning  if  the  time  step  size  must  be  decreased 
dramatically  to  meet  the  error  criterion.  But  in  the  gun  code  there  are 
several  discontinuities.  For  example,  when  the  chamber  pressure  first 
exceeds  the  shot  start  pressure,  the  ordinary  differential  equations  for 
the  projectile  motion  are  suddenly  changed.  The  code  normally  has  to 
cut  back  on  the  time  step  to  integrate  through  this  and  similar  changes. 
The  normal  warning  messages  are  suppressed  unless  the  write  option  is 
turned  on  (KWRITE-1). 

The  routine  TSTEP  actually  takes  the  integration  steps.  A  number 
of  tests  have  been  implemented  here.  For  example,  Eq.  [4]  for  the 
volume  of  the  liquid  reservoir  becomes  singular  as  V^  approaches  zero. 

So  when  V^  becomes  less  then  the  initial  volume  divided  by  a  thousand, 
the  region  is  closed.  That  is,  region  1  is  ignored  in  the  integration 
from  then  on.  The  error  resulting  from  ignoring  the  small  amount  of 
remaining  propellant  is  negligable.  If  becomes  negative,  the  time 
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step  was  too  large,  and  it  must  be  reduced  and  the  time  step  taken  over. 
There  is  a  similar  problem  with  region  4,  since  the  projectile  may  start 
at  the  throat,  and  V4  will  be  zero.  In  this  case  region  4  starts  out  as 
closed.  When  the  projectile  is  at  least  .1  cm.  down  the  gun  tube,  the 
region  is  opened.  Before  this,  region  4  is  considered  as  an  extension 
of  region  3.  Also,  diagnostic  printouts  of  the  quantities  of  interest 
may  be  made  after  each  time  step  (KWRITE-1) . 

There  are  two  main  integration  routines  included  In  EPISODE.  The 
first  is  the  Adams  method,  suitable  for  non-stiff  problems  (METH-1) . 

The  second  is  the  backward  differentiation  method,  suitable  for  stiff 
problems  (METH-2).  We  have  normally  used  METH-2.  However,  the 
governing  equations  are  not  very  stiff,  and  METH-1  can  also  be  used. 

The  run  times  are  about  equivalent  for  the  two  methods. 

There  are  also  several  methods  for  the  nonlinear  system  solution 
required  at  each  step  of  the  integration.  We  have  used  the  Newton 
method  with  an  internally  computed  finite  difference  approximation  to 
the  Jacobian.  So  the  only  subroutine  that  needs  to  be  supplied  by  the 
user  is  DIFFUN.  Given  a  set  of  values  for  the  unknowns,  DIFFUN  computes 
the  time  derivatives  of  the  unknowns. 

Subroutine  DIFFUN  in  this  case  only  calls  one  of  two  possible 
subroutines  (FDR0P1  or  FDR0P2).  These  routines  correspond  to  the 
instantaneous  burning  option  and  the  droplet  burning  option.  Since  the 
governing  equations  are  different,  It  is  easier  to  set  these  up  as 
separate  subroutines. 

The  main  routine  RLGD  reads  in  the  input  data.  Initial  values  are 
assigned  in  INITAL.  The  subroutine  INTEG  actually  calls  the  integrator 
DRIVE.  If  a  repeat  option  has  been  specified,  INTEG  will  reset  the 
initial  conditions  and  repeat  until  convergence  has  occured. 

Most  of  the  subroutines  are  for  the  different  options.  To  make  it 
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easier  to  add  or  change  options,  some  of  the  input  is  in  separate 
subroutines.  For  instance,  subroutines  VENT1 ,  VENT2,  and  VENT3  read  in 
the  required  data  for  the  three  vent  options.  These  routines  also 
compute  the  vent  area  for  a  given  piston  travel.  So  if  a  new  vent 
option  is  desired,  it  is  only  necessary  to  add  a  new  subroutine  and 
change  the  few  lines  that  call  the  VENT  subroutines.  The  subroutines 
for  the  other  options  are  set  up  similarly. 

Subroutines  CAPTION  and  OUT  control  the  output  files.  Subroutine 
OUTGRA  creates  a  graphics  file,  so  that  any  of  the  quantities  computed 
may  later  be  graphed. 


XXIII.  BRL  3<"MM  GUN  FIXTURE 

A  set  of  experimental  measurements  have  been  made  on  a  30mm 
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regenerative  liquid  propellant  gun  at  the  BRL.  ’  The  liquid 
reservoir  pressure,  combustion  chamber  pressure,  piston  travel,  and 
projectile  velocity  have  been  measured.  The  data  has  been  filtered  to 
remove  the  acoustic  oscillations.  A  number  of  cases  have  been  measured 
for  a  2/3  charge  and  a  1/3  charge.  The  gun  has  not  been  fired  with  a 
full  charge.  The  injection  thickness  of  the  annular  jet  has  -not  been 
optimized,  and  the  performance  of  the  gun  is  low.  In  this  paper  only 
one  experiment  will  be  considered  (Round  8,  2/3  charge). 

Figure  3  shows  a  scale  model  drawing  of  the  liquid  reservoir 

before  firing.  The  piston  is  flush  against  a  crash  ring  in  the 

combustion  chamber  (not  shown) .  The  crash  ring  is  for  protection  in 

case  of  a  piston  reversal.  An  0-ring  initially  seals  the  liquid 

reservoir.  The  block  should  be  exactly  at  the  end  of  the  back  taper  of 

the  central  bolt.  However,  when  the  fixture  is  assembled,  the  block 

2  3 

actually  protudes  about  an  eighth  of  the  inch  forward.  When  the 
liquid  propellant  is  loaded,  the  reservoir  is  presurized  to  7  MPa.  This 
is  enough  to  push  the  block  back,  so  the  block  will  be  only  slightly 
forward  of  the  back  taper. 
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30mm  Gun  -  2/3  charge 


Figure  3.  Scale  Model  Drawing  of  the  Liquid  Reservoir,  Initial  Position 
(after  Pre-pressurization) . 

Figure  4  shows  the  experimental  pressure  profiles.  The  primer 
causes  the  initial  pressure  rise  in  the  chamber.  The  chamber  pressure 
then  remains  constant  while  the  liquid  pressure  slowly  rises.  Both  the 
piston  and  the  block  are  in  motion.  When  the  Belleville  springs  bottom 
out,  the  momentum  of  the  piston  causes  a  pressure  rise  in  the  reservoir. 
This  sets  up  oscillations  in  the  liquid  pressure.  The  injected  liquid 
begins  to  combust  and  the  chamber  and  liquid  pressures  increase  rapidly. 

Figure  5  shows  the  experimental  piston  travel.  The  piston  moves 
relatively  slowly  at  first.  When  the  Belleville  springs  bottom  out,  the 
pressure  rise  in  the  liquid  causes  the  piston  to  stop.  After  this  the 
main  injection  stroke  takes  place.  Figure  6  gives  the  experimental 
projectile  velocity.  This  is  only  accurate  over  a  period  of  time  near 
the  initial  projectile  motion.  Tue  muzzle  velocity  (measured 
separately)  is  around  1020  m/s. 
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Figure  6.  Experimental  Projectile  Velocity. 


*1  ^ 

In  a  previous  report  the  data  was  analyzed  using  an  inverse 
code.  Values  of  the  discharge  coefficient  and  the  liquid  accumulation 
were  derived.  The  discharge  coefficient  started  small  and  took  several 
milliseconds  to  rise  to  a  value  near  one  (see  Figure  7).  It  is 
sometimes  more  convenient  to  consider  the  discharge  coefficient  versus 
the  piston  travel  (Figure  8).  At  early  times  the  derived  discharge 
coefficient  is  not  expected  to  be  accurate.  A  small  error  in  the  piston 
position  can  lead  to  a  relatively  large  change  in  the  vent  area.  It  is 
interesting  that  even  after  the  piston  is  past  the  front  taper  (piston 
travel  -  1.32  cm.)  the  discharge  coefficient  is  still  small. 

More  recently**  an  attempt  was  made  to  duplicate  this  behavior 
using  transient  models  of  the  liquid  injection  (lumped  parameter,  one¬ 
dimensional,  and  two-dimensional).  The  two-dimensional  model  has  only 
been  applied  to  simplified  cases  where  the  central  bolt  does  not  have  a 
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taper.  All  the  models  indicated  a  rapid  rise  to  a  discharge  coefficient 
near  one.  This  led  to  a  more  careful  look  at  the  inverse  procedure. 

The  derivation  of  the  discharge  coefficients  depends  on  very  accurate 
measurements  of  the  piston  travel.  In  fact,  the  discharge  coefficients 
depend  on  the  piston  velocity,  so  the  data  must  be  numerically 
differentiated.  However,  the  block  motion  has  not  been  measured,  so 
this  must  be  estimated.  And  the  recoil  of  the  gun  has  not  been 
measured.  Since  only  the  relative  piston  motion  is  measured,  and  this 
must  be  scaled  against  the  expected  maximum  piston  travel,  an  error 
anywhere  during  the  gun  firing  will  effect  the  entire  profile.  Also,  it 
is  known  that  the  mini- transducer  used  in  the  liquid  reservoir  is  less 
accurate  than  the  transducers  in  the  combustion  chamber.  If  the  flow 
was  actually  steady  state,  the  liquid  pressure  would  just  be  the  chamber 
pressure  times  the  hydraulic  difference.  However,  the  oscillations  in 
the  liquid  pressure  lead  to  substantial  differences  in  the  ratio  of  the 
liquid  to  the  chamber  pressure,  and  the  accuracy  of  the  liquid  pressure 
measurement  is  unknown.  Therefore  the  derivation  of  the  discharge 
coefficient  is  here  given  less  weight,  and  an  attempt  is  made  to  model 
the  experiment  assuming  a  constant  large  discharge  coefficient. 

Appendix  B  shows  the  input  for  the  code  for  this  problem.  The 
numbers  and  labels  at  the  left  are  read  in  by  the  code.  The  comments  at 
the  right  are  for  Identification  by  the  user  and  do  not  effect  the 
actual  code.  The  corresponding  output  follows. 

The  first  line  is  merely  a  label  for  the  problem.  It  lists  the 
filename  of  the  input  job  stream  and  a  brief  description  of  the  problem. 

The  initial  offset  of  the  projectile  and  the  total  distance 
traveled  by  the  projectile  before  muzzle  exit  are  given.  The  diameter 
of  the  gun  tube  is  given.  The  projectile  and  piston  weights  are  entered 
(measured) . 

The  initial  volumes  of  the  reservoir  and  chamber  are  given.  The 


50 


reservoir  volume  is  computed  from  the  drawings.  The  crash  ring  in  the 
chamber  makes  computing  the  volume  difficult,  so  the  initial  volume  is 
measured.  The  initial  gun  tube  volume  will  be  computed  from  the 
offset  and  the  gun  tube  diameter. 

The  initial  areas  of  the  reservoir  and  chamber  are  given  (from  the 
engineering  drawings) .  The  gun  tube  area  is  computed  from  the  gun  tube 
diameter. 

The  vent  option  is  chosen.  In  this  case  the  actual  piston  option 
is  picked  (VENT3).  This  requires  a  table  of  the  bolt  radii  versus 
distance  (from  drawings).  The  zero  distance  is  che  initial  position  of 
the  front  of  the  piston.  The  vent  area  is  the  area  between  the  front  of 
the  piston  and  the  bolt.  The  area  of  this  central  hole  in  the  piston  is 
given,  as  well  as  the  area  of  the  grease  dyke.  The  O-ring  originally  is 
between  0.127  and  0.363  cm.  measured  from  the  initial  piston  position. 

As  the  0-ring  is  uncovered,  the  high  liquid  pressure  will  push  out  the 
0-ring  before  the  piston  has  moved  out  to  0.363  cm.  The  exact  behavior 
of  the  vent  is  not  known  here.  As  an  approximation,  the  bolt  radius  is 
assumed  to  vary  linearly  from  the  beginning  to  the  end  of  the  initial 
location  of  the  0-ring. 

The  Belleville  spring  option  is  chosen.  The  force  exerted  by  the 
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Belleville  springs  has  been  measured  as  a  function  of  distance.  The 
springs  bottom  out  at  0.422  cm.  The  pre-presurization  is  sufficient  to 
compress  the  springs  0.184  cm.,  so  the  block  only  moves  0.238  cm.  during 
the  gun  firing.  The  surface  area  of  the  block  (from  the  drawings)  and 
the  weight  of  the  block  (measured)  must  also  be  known. 

For  lack  of  better  information,  the  piston  resistance  is  set  equal 
to  zero.  The  discharge  coefficient  into  the  chamber  is  set  at  a 
constant  0.95.  This  agrees  well  with  higher  dimensional  transient 
models.^  The  discharge  coefficient  into  the  gun  tube  is  chosen  as  one 
(no  losses)  . 
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The  flow  into  the  chamber  is  modeled  as  steady  state  Bernoulli 
flow  (FLUX1) .  There  is  only  one  vent  hole.  The  piston  thickness  is 
irrelevant  for  this  option.  The  flow  into  the  gun  tube  is  steady  state 
isentropic  flow  (FLUX1) . 

The  shot  start  pressure  is  taken  as  50  MPa.  Previously  a  value  of 

68  MPa  has  been  used.  An  analysis  of  the  radar  traces  of  the  projectile 

23 

indicate  the  projectile  moves  somewhat  sooner.  The  initial  projectile 
motion  may  not  be  reproducible,  and  there  is  some  evidence  that  the 

projectile  sometimes  moves  a  few  centimeters,  stops,  and  then  later 

0  / 

moves  again.  After  shot  start  the  resistance  pressure  is  taken  as  5 
MPa.  This  is  a  typical  value  from  solid  propellant  gun  modeling. 

Next  the  physical  properties  of  the  propellant  HAN1846  are  given. 

g 

The  density  at  atmospheric  pressure  has  been  measured  and  the  bulk 
modulus  has  been  fitted. ^  The  chemical  energy,  the  ratios  of  specific 
heats,  the  molecular  weight  of  the  propellant  final  produces  gas,  and 

Q 

the  covolume  have  been  computed  using  BLAKE.  The  surface  tension  and 
the  kinematic  viscosity  of  the  liquid  are  not  actually  required,  but 
are  used  to  compute  the  Weber  number,  the  Reynold's  number,  and  the 
third  Lagrange  number. 

The  liquid  is  pre -presurized  to  7.0  MPa.  The  initial  chamber 
pressure  is  one  atmosphere. 

As  a  first  approximation,  the  liquid  is  assumed  to  combust 
instantaneously  as  it  enters  the  combustion  chamber  (DR0P1).  Later  the 
effects  of  accumulation  will  be  studied. 

The  primer  is  assumed  to  be  injected  in  the  form  of  hot  gas 
(PRIM3).  From  the  experimental  chamber  pressure,  once  the  pressure 
starts  to  rise,  it  takes  about  2.5  ms.  to  reach  a  pressure  of  13  MPa. 

The  chamber  pressure  then  stays  almost  constant  for  another  1.5  ms.  The 
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initial  pressure  rise  is  assumed  to  be  due  solely  to  the  primer.  So  the 
injection  time  is  taken  to  be  2.5  ms.  The  primer  mass  is  chosen  as  1.5 
grams  in  order  to  obtain  the  proper  pressure.  The  actual  primer  is  3.0 
grams.  But  since  the  primer  is  injected  through  a  relatively  long 
narrow  tube,  there  should  be  large  energy  losses. ^  The  code  does  not 
model  these  losses  directly,  so  a  smaller  amount  of  primer  must  be 
chosen. 

The  heat  loss  to  the  combustion  chamber  is  ignored.  Heat  loss  to 
the  gun  tube  walls  is  calculated.  The  gun  tube  wall  temperature  is 
taken  to  be  300X.  The  heat  loss  algorithm  in  the  code  is  not  adjusted. 

The  air  shock  in  front  of  the  projectile  is  calculated.  The  air 
in  the  tube  is  initially  assumed  to  be  at  1  atm.  and  300K.  The  values 
of  the  molecular  weight  and  ratio  of  specific  heats  are  taken  at  these 
conditions . 

Since  this  is  a  low  performance  gun,  the  choice  of  the  pressure 
distribution  in  the  gun  tube  is  not  critical.  The  modified  Lagrange 
distribution  is  chosen  (TUBE2) .  This  does  take  into  account  the 
velocity  at  the  entrance  to  the  tube,  but  is  still  easy  to  calculate. 
There  is  no  water  buffer  in  this  fixture  (BUFFI). 

The  code  will  print  out  results  every  0.1  ms.  (TINC).  Because  the 
code  must  often  change  the  time  step,  it  is  more  efficient  to  restrict 
the  maximum  rime  step  (HTOF;.  The  error  controls  EPS  and  SREC  are  given 
typical  values . 

The  integration  method  flag  MF  is  set  to  22  (backwards 
differentiation  formulas  with  an  internally  computed  Jacobian) .  KWRITE 
is  set  to  zero  to  eliminate  diagnostic  messages.  A  time  limit  TMAX  is 
set.  If  the  code  takes  longer  than  TMAX  seconds  to  execute,  the  code 
will  stop  and  write  the  usual  summary  pages  so  all  the  information  will 
not  be  lost. 
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The  code  is  only  to  be  integrated  once  (REP1)  and  the  chamber 
pressure  will  be  computed  normally  (CHAM1). 

Appendix  B  also  gives  the  output  for  this  problem.  The  code  was 
run  on  the  BRL  Cray-2.  The  input  is  echoed  back,  along  with  other 
calculated  initial  values.  Once  the  integration  begins,  output  is  given 
every  0.1  ms.  The  pressures,  piston  travel  and  velocity,  and  projectile 
travel  and  velocity  are  given.  At  the  end  of  the  integration,  a  summary 
page  is  given.  A  second  page  shows  when  various  maximum  values  are 
achieved.  There  are  many  other  variables  of  interest  in  the  simulation. 
As  the  integration  proceeds,  these  are  written  onto  other  files,  and  can 
be  appended  to  the  output.  In  Appendix  B,  only  the  minimal  output  is 
given. 

Figure  9  shows  the  resulting  chamber  pressure  profile.  Because  of 
the  instantaneous  burning  option,  the  pressure  starts  to  rise  as  soon  as 
the  vent  opens.  The  only  way  to  delay  the  pressure  rise  would  be  to 
delay  the  opening  of  the  vent. 

To  make  comparisons  more  easily,  the  model  results  are  shifted 
(Figures  10-13).  The  curves  are  moved  over  so  that  the  model  chamber 
pressure  reaches  the  shot  start  pressure  (50  MPa)  at  the  same  time  as 
the  experimental  chamber  pressure.  All  the  model  curves  are  moved  the 
same  amount.  After  the  initial  disturbances  caused  by  the  Belleville 
springs,  the  agreement  in  the  pressure  curves  is  reasonable.  The 
agreement  in  piston  travel  is  poor.  The  model  projectile  velocities  are 
consistently  high,  due  to  the  higher  chamber  pressure,  and  the  predicted 
muzzle  velocity  is  about  6%  high. 

A  droplet  option  was  also  calculated  to  see  if  improved  agreement 
could  be  obtained.  The  inverse  code  was  not  used  to  obtain  liquid 
accumulation  because  of  doubts  concerning  the  accuracy  of  the  piston 
travel.  Instead,  the  drop  size  versus  piston  travel  table  was  adjusted 
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Figure  10.  Chamber  Pressure  -  Centered.  Experiment  (line) 
Model  -  Instantaneous  Burning  -  Cn  -  0.95  (dot). 


Figure  13.  Projectile  Velocity  -  Centered.  Experiment  (line). 

Model  -  Instantaneous  Burning  -  Cjj  -  0.95  (dot). 

by  trial  and  error  until  reasonable  agreement  in  the  chamber  pressures 
was  obtained.  The  DR0P1  option  then  is  replaced  by  the  DR0P3  option. 
This  is  the  only  change  made  in  the  code .  The  input  and  the  summary 
page  of  the  output  are  given  in  Appendix  C.  The  results  (Figures  14-17) 
no  longer  require  centering,  since  the  delay  in  the  experimental 
pressure  rise  is  now  being  modeled. 

The  agreement  in  the  combustion  chamber  pressures  is  much  better. 
The  oscillations  in  the  liquid  pressure  are  now  reproduced.  However,  in 
the  model  the  Belleville  springs  bottom  out  too  early.  This  can  be  seen 
in  the  liquid  pressure  profile,  where  the  pressure  bump  comes  too  early, 
and  in  the  piston  travel  profile,  where  the  piston  stops  sooner.  The 
later  agreement  in  the  piston  travel  may  or  may  not  be  coincidental. 

The  muzzle  velocity  is  still  3%  high. 
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Liquid  pressure  (MPa) 


Figure  15.  Liquid  Pressure.  Experiment  (line) 
Model  -  Droplet  Formation  -  Cn  -  0.95  (dot). 


The  experiment  indicates  that  the  block  should  move  further  before 
the  Belleville  springs  bottom  out.  So  calculations  were  made  where  the 
block  is  assumed  to  move  2  mm.  further.  This  would  indicate  a  small 
error  in  the  measurement  of  the  spring  behavior.  Note  that  the 
experimental  piston  travel  must  now  be  scaled  to  a  slightly  longer  total 
piston  travel. 

The  model  with  instantaneous  burning  is  calculated  first.  The 
input  is  in  Appendix  D,  and  Figures  18-21  show  the  results.  The 
agreement  is  now  very  poor.  The  piston  moves  another  2  mm.  before  the 
block  stops  moving.  Therefore  the  vent  area  is  larger,  and  the 
injection  starts  more  rapidly.  Because  of  the  assumption  of 
instantaneous  burning,  the  initial  small  difference  in  the  injection 
rate  grows  as  time  goes  on.  The  unfortunate  fact  here  is  that  the 
initial  piston  position  must  be  known  very  accurately  to  model  this  gun 
fixture  properly.  Fortunately  later  gur  fixtures  do  not  use  Belleville 
springs  or  an  O-ring,  but  a  metal- to-metal  seal,  and  the  vent  area 
versus  piston  travel  is  much  better  known. 

Next  droplet  burning  was  assumed,  and  again  the  droplet  size  table 
was  chosen  so  that  the  chamber  pressure  profile  was  reasonably  accurate 
(Appendix  E,  Figures  22-25).  The  agreement  in  the  piston  travel  is 
still  rather  poor.  Despite  experimentation  with  the  code,  it  proved 
impossible  to  get  reasonable  agreement  with  the  piston  travel  under  the 
assumption  that  the  discharge  coefficient  was  uniformly  equal  to  0.95. 

So  an  attempt  was  made  to  get  agreement  in  the  piston  travel 
making  minimal  changes  in  the  discharge  coefficient.  This  proved  to  be 
possible  by  decreasing  the  discharge  coefficient  when  the  piston  was 
over  the  front  taper  of  the  bolt  and  increasing  it  back  to  0.95  when  the 
piston  moved  to  the  flat  part  of  the  bolt  (Appendix  F,  Figures  26-29). 
The  agreemert  with  the  chamber  pressure  and  the  piston  travel  is  now 
very  good. 
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Figure  18.  Chamber  Pressure  -  Centered.  Experiment  (line) 

Model  -  Instantaneous  Burning  -  CQ  -  0.95  -  Long  Belleville  (dot). 


Figure  19.  Liquid  Pressure  *  Centered.  Experiment  (line). 

Model  -  Instantaneous  Burning  -  CD  -  0.95  -  Long  Belleville  (dot). 
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Figure  20.  Piston  Travel  -  Centered.  Experiment  (line). 

Model  -  Instantaneous  Burning  -  Cp  -  0.95  -  Long  Belleville  (dot). 


Figure  21.  Projectile  Velocity  -  Centered.  Experiment  (line). 
Model  -  Instantaneous  Burning  -  C^  —  0.95  -  Long  Belleville  (dot) 
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Figure  22.  Chamber  Pressure.  Experiment  (line). 

Model  -  Droplet  Formation  -  -  0.95  -  Long  Belleville  (dot) 


Figure  23.  Liquid  Pressure.  Experiment  (line). 

Model  -  Droplet  Formation  -  Cn  -  0.95  -  Long  Belleville  (dot) 
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Figure  24.  Piston  Travel.  Experiment  (line). 

Model  *  Droplet  Formation  -  Cq  -  0.95  -  Long  Belleville  (dot). 


Figure  25.  Projectile  Velocity.  Experiment  (line). 

Model  -  Droplet  Formation  -  Cq  -  0.95  -  Long  Belleville  (dot). 
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Figure  26.  Chamber  Pressure.  Experiment  (line). 

Modol  Droplet  Formation  -  Cq  Adjusted  -  Long  Bellevil 


CD 

C. 

o 

cn 

cn  150.0 

® 

c 

CL 


Figure  27.  Liquid  Pressure, 
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Figure  28.  Piston  Travel.  Experiment  (line). 

Model  -  Droplet  Formation  -  CD  Adjusted  -  Long  Belleville  (dot 


Figure  29.  Projectile  Velocity.  Experiment  (line). 

Model  -  Droplet  Formation  -  CD  Adjusted  -  Long  Belleville  (dot) 


Figure  30.  Discharge  Coefficient.  Derived  from  Experiment  (line). 
Model  -  Droplet  Formation  -  Cq  Adjusted  -  Long  Belleville  (dot) . 


Figure  30  shows  a  comparison  of  the  present  profile  for  the 
discharge  coefficient  with  the  experimentally  derxved  values.  The 
differences  are  due  partly  to  the  small  differences  in  the  piston 
travel,  since  the  discharge  coefficient  depends  on  the  slope  of  the 
piston  travel.  The  early  differences  are  due  to  the  differences  in  the 
liquid  pressure.  Since  the  liquid  pressure  transducer  is  known  to  be 
less  accurate,  and  since  the  ratio  of  the  liquid  to  chamber  pressure  is 
equal  to  the  hydraulic  difference  only  at  steady  state,  the  model  liquid 
pressures  are  still  plausible. 

To  summarize,  the  last  model  shows  very  good  agreement  with 
experiment.  The  assumption  must  be  made  that  the  discharge  coefficient 
is  small  when  the  piston  is  over  the  front  taper  of  the  bolt.  This 
occurs  right  after  the  Belleville  springs  bottom  out.  The  liquid 
pressure  oscillations  could  lead  to  turbulence  and/or  flow  separation  at 
the  exit.  At  this  point,  the  opening  is  essentially  a 
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converging/diverging  nozzle.  The  earlier  lumped  parameter  and  one¬ 
dimensional  transient  models  would  not  be  able  to  predict  this  type  of 
behavior,  and  the  two-dimensional  model  has  not  been  applied  to  a 
problem  with  a  tapered  bolt. 

The  model  does  not  reproduce  the  pressure  behavior  after  the 
maximum  pressure.  At  this  stage  the  piston  is  near  the  end  of  its 
stroke.  Hydrodynamic  effects  may  be  important.  Since  the  model 
pressure  is  uniformly  higher  than  the  experimental  pressure  after  the 
maximum,  there  are  probably  loss  terms  not  implemented  in  the  model. 
The  muzzle  velocity  is  also  slightly  too  large. 

XXIV.  CONCLUSIONS 


An  update  is  given  to  a  lumped  parameter  model  for  regenerative 
liquid  propellant  guns.  While  the  basic  structure  is  the  same,  many 
more  options  have  been  added.  A  complete  discussion  of  the  governing 
equations  and  the  different  options  is  given.  The  code  has  been  used  to 
test  the  effect  of  various  assumptions,  so  additional  options  will  be 
added  as  needed. 

The  model  results  are  compared  with  experimental  data  from  a  30mm 
gun  fixture.  A  discussion  is  given  of  the  procedure  for  choosing  a  set 
of  input  parameters.  The  code  predicts  the  overall  behavior  of  the  gun 
fixture  extremely  well.  Some  of  the  details  are  not  reproduced,  and  the 
probable  shortcomings  in  the  code  are  discussed. 
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GLOSSARY 


o 

Area  of  the  liquid  reservoir,  cm  . 

2 

A^  Area  of  the  chamber,  cm  . 

2 

A^  Area  of  the  gun  tube ,  cm  . 

2 

A^£  Area  of  the  piston  in  the  buffer,  cm  . 

2 

Aj^  Area  of  the  block  in  the  reservoir,  cm  . 

2 

Ag  Area  of  the  grease  dyke  around  an  annular  piston,  cm  . 

2 

A-^  Area  of  the  central  hole  in  the  annular  piston,  cm  . 

2 

Av  Area  of  the  piston  vents,  cm  . 

b  Covolume,  cm  /g, 

c^  Speed  of  sound  in  the  liquid  in  the  reservoir,  cm/s. 

C3  Speed  of  sound  in  the  mixture  in  the  chamber,  cm/s . 

c^  Speed  of  sound  in  the  mixture  in  the  gun  tube,  cm/s. 

CL3  Speed  of  sound  in  the  liquid  in  the  chamber,  cm/s . 

Speed  of  sound  in  the  liquid  in  the  gun  tube,  cm/s. 
c Speed  of  sound  in  the  gas  in  the  chamber,  cm/s. 
c Average  speed  of  sound  in  the  gas  in  the  gun  tube,  cm/s. 
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c^£  Speed  of  sound  in  the  buffer  liquid,  cm/s. 

Cp  Specific  heat  at  constant  pressure,  j/g-K. 

cv  Specific  heat  at  constant  volume,  j/g-K. 

Cp  Discharge  coefficient  for  the  mass  flux  through  the  piston. 

Cp'  Discharge  coefficient  for  the  mass  flux  into  the  gun  tube, 

d^  Diameter  of  the  gun  tube,  cm. 

e-^  Chemical  energy  of  the  liquid,  j/g. 

Internal  energy  of  the  gas  in  the  chamber,  j/g. 
e^  Internal  energy  of  the  gas  in  the  gun  tube,  j/g. 

Ejj^  Total  internal  energy  in  the  liquid  reservoir,  j. 

EL3  Total  internal  energy  in  the  liquid  in  the  chamber,  j. 

E^4  Total  internal  energy  in  the  liquid  in  the  gun  tube,  j. 

Eq^  Total  internal  energy  in  the  gas  in  the  chamber,  j. 

Eq4  Total  internal  energy  in  the  gas  in  the  gun  tube,  j. 

Ep  Total  internal  energy  in  the  unburned  primer,  j. 

EKL4  Kinetic  energy  of  the  liquid  in  the  gun  tube,  j. 

EKq^  Kinetic  energy  of  the  gas  in  the  gun  tube,  j. 
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EKpj 

EKPS 

EH3 

eh4 


EFPJ 


EF 


ps 


Lbk 


’Ll 


*L3 


L4 


*G3 


G4 


K 

K 


Kinetic  energy  of  the  projectile,  j. 

Kinetic  energy  of  the  piston,  j. 

Total  heat  loss  to  the  chamber  walls,  j. 

Total  heat  loss  to  the  gun  tube  walls,  j. 

Total  energy  loss  due  to  projectile  friction,  j . 

Total  energy  loss  due  to  piston  friction,  j . 

Total  energy  o f  the  system,  j. 

2 

Force  exerted  by  the  Belleville  springs,  MPa-cm  . 

7  2 

Conversion  constant  -  10  g/MPa-cm-s  . 

Liquid  enthalpy  in  the  reservoir,  j/g. 

Liquid  enthalpy  in  the  combustion  chamber,  j/g- 

Liquid  enthalpy  in  the  gun  tube,  j/g. 

Gas  enthalpy  in  the  chamber,  j/g. 

Gas  enthalpy  in  the  gun  tube,  j/g. 

2 

Heat  transfer  coefficient  to  the  chamber  walls,  j/cm  -K-s. 

2 

Heat  transfer  coefficient  to  the  gun  tube  walls,  j/cm  -K-s. 
Bulk  modulus,  MPa. 
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K1 

k2 

L 


m13 

m34 

"^3 

®bf 

m4 

mL3 

mL4 

mG3 

mG4 

m3 

M4 

ML3 


Bulk  modulus  at  zero  pressure,  MPa. 

Derivative  of  the  bulk  modulus . 

Piston  thickness,  cm. 

Mass  flux  from  the  reservoir  into  the  chamber,  g/s. 

Mass  flux  from  the  chamber  into  the  gun  tube,  g/s. 

Mass  flux  of  the  primer  into  the  chamber,  g/s. 

Mass  flux  out  of  the  buffer,  g/s. 

Chamber  rate  of  production  of  gas  by  combustion,  g/s. 

Gun  tube  rate  of  production  of  gas  by  combustion,  g/s. 
Total  rate  of  production  of  liquid  in  the  chamber,  g/s. 

Total  rate  of  production  of  liquid  in  the  gun  tube,  g/s. 

Total  rate  of  production  of  gas  in  the  chamber,  g/s. 

Total  rate  of  production  of  gas  in  the  gun  tube,  g/s. 

Total  mass  in  the  reservoir,  g. 

Total  mass  in  the  chamber,  g. 

Total  mass  in  the  gun  tube,  g. 

Liquid  mass  in  the  chamber,  g. 
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M 


L4 


M, 


G3 


M, 


G4 


*P 

MJ 


Mbk 

Kps 

Mg 

Pi 


P4 

PL 

PR 

Pt 

PPj 


ps 


Liquid  mass  in  the  gun  tube,  g. 

Gas  mass  in  the  chamber,  g. 

Gas  mass  in  the  gun  tube,  g. 

Unburned  primer  mass,  g. 

Total  propellant  mass  in  the  system,  g. 
Mass  of  the  block,  g. 

Mass  of  the  projectile,  g. 

Mass  of  the  piston,  g. 

Molecular  weight  of  the  gas,  g/mole. 
Pressure  in  the  liquid  reservoir,  MPa. 
Pressure  in  the  chamber,  MPa. 


Space  mean  pressure  in  the  gun  tube,  MPa. 
Pressure  at  the  gun  throat,  MPa. 

Pressure  at  the  base  of  the  projectile,  MPa. 
Pressure  at  the  gun  tube  entrance  -  p^,  MPa. 
Projectile  resistance  pressure,  MPa. 

Piston  resistance  pressure,  MPa. 
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Ps 

Pbf 

Pr 

Qc 

Qw 

Re 

Rs 

*u 


s 


s 


s 


bk 

Pj 

ps 


S 

S 


3 

4 


t 


T 

T 

T 


3 

4 

o 


Air  shock  pressure  ahead  of  the  projectile,  MPa. 
Pressure  in  the  buffer,  MPa. 


Prandtl  number. 

3 

Heat  loss  to  the  chamber  walls,  j/cm  -s. 

3 

Heat  loss  to  the  gun  tube  walls,  j/cm  -s. 
Reynold's  number. 

Specific  gas  constant,  j/g-K. 

Universal  gas  constant  -  8.318  j/mole-K. 
Block  travel,  cm. 

Projectile  travel,  cm. 

Piston  travel,  cm. 


Chamber  mass  production  rate  times  enthalpy  change  ,  j/s. 
Gun  tube  mass  production  rate  times  enthalpy  change  ,  j/s. 

Time,  s. 


Temperature  in  the  combustion  chamber,  K. 
Average  temperature  in  the  gun  tube,  K. 
Isoc'  >ric  temperature  of  the  propellant,  K. 
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Tt  Temperature  at  the  gun  tube  entrance,  K. 

Tc  Temperature  of  the  chamber  walls,  K. 

Tw  Temperature  of  the  gun  tube  walls,  K. 

Injection  velocity  of  the  liquid  into  the  chamber,  cm/s. 
Velocity  of  the  block,  cm/s. 

Vpj  Velocity  of  the  projectile,  cm/s. 

v  Velocity  of  the  piston,  cm/s. 

3 

V^  Volume  of  the  liquid  reservoir,  cm  . 

O 

V3  Volume  of  the  chamber,  cm  . 

V^  Volume  of  the  gun  tube  behind  the  projectile,  cm  . 

VL3  Volume  of  the  liquid  in  the  chamber,  cm  . 

3 

Vla  Volume  of  the  liquid  in  the  gun  tube,  cm  . 

Vq^  Volume  of  the  gas  in  the  chamber,  cm  . 

3 

^G4  Volu  le  of  the  gas  in  the  gun  tube,  cm  . 

3 

V^f  Volume  of  the  buffer,  cm  . 

vt  Fluid  velocity  at  the  gun  tube  entrance,  cm/s. 

Xg  Distance  from  the  tube  entrance  to  the  projectile  base,  cm. 
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7 


Ratio  of  specific  heats. 


£3  Porosity  of  the  mixture  in  the  combustion  chamber. 

£4  Porosity  of  the  mixture  in  the  gun  tube. 

X  Impetus  of  the  propellant,  j/g. 

p  Viscosity  of  the  gas,  poise  -  g/cm-s. 


Po 
P 1 

P3 

PU 

PL3 
P 1A 

PC  3 


Viscosity  of  the  gas  at  the  wall  temperature,  poise. 

3 

Liquid  density  at  zero  pressure,  g/cm  . 

3 

Liquid  density  in  the  reservoir,  g/cm  . 

3 

Mixture  density  in  the  chamber,  g/cm  . 

3 

Mixture  density  in  the  gun  tube,  g/cm  . 

3 

Liquid  density  in  the  chamber,  g/cm  . 

Liquid  density  in  the  gun  tube,  g/cm  . 

3 

Gas  density  in  the  chamber,  g/cm  . 


PGU 

pbf 


3 

Gas  density  in  the  gun  tube,  g/cm  . 

3 

Liquid  density  in  the  buffer,  g/cm  . 
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APPENDIX  A 


The  computer  code  is  long  (around  100  pages) ,  so  a  complete 
listing  is  not  given.  Instead,  only  the  description  of  the 
input  options  from  the  beginning  of  the  code  is  given.  A 
complete  listing  of  the  code  (tape,  disk,  or  hardcopy)  may 
be  obtained  from  the  author. 


Q  A******************************************** 

C  FEBRUARY  23,  1988. 

C  DESCRIPTION  OF  INPUT 

C 

C  FILE  1:  PROBLEM  TITLE 

C  80  ALPHANUMERIC  CHARACTERS. 

C  WILL  BE  TITLE  ON  GRAPHICS  FILES. 

C 

C  FILE  2:  OFFSET,  TRAVEL 

C  OFFSET=DISTANCE  FROM  END  OF  TUBE  TO  PROJECTILE  (CM) . 

C  TRAVEL=TOTAL  PROJECTILE  TRAVEL  (CM) . 

C 

C  FILE  3 : D4 

C  D4  =  GUN  TUBE  DIAMETER  (CM) . 

C 

C  FILE  4:  PJWT 

C  PJWT  =  PROJECTILE  WEIGHT  (GM) . 

C 

C  FILE  5:  PSWT 

C  PSWT  =  PISTON  WEIGHT  (GM) . 

C 

C  FILE  6:  VII, V3I 

C  VII  =  INITIAL  LIQUID  RESERVOIR  VOLUME  (CM**3). 

C  V3I  =  INITIAL  COMBUSTION  CHAMBER  VOLUME  (CM**3). 

C 

C  FILE  7:  A1 , A3 

C  A1  =  AREA  OF  LIQUID  SIDE  PISTON  FACE  (INCLUDING  VENTS)  (CM**2) . 

C  A3  =  AREA  OF  CHAMBER  SIDE  PISTON  FACE  (INCLUDING  VENTS)  (CM**2). 

C 

C  FILE  8:  TVENT  -  VENT  OPTION  SWITCH 

C 

C  IF  VENT1  THEN  PISTON  WITH  HOLES. 

C  AREA  OF  HOLES  MAY  VARY  WITH  PISTON  TRAVEL. 

C  ENTER  VENT  AREA  TABLE. 

C  FILE  8.1:  NVENT 

C  NVENT  =  NO.  OF  TABLE  ENl..xES. 

C  FILE  8.2...: SVT ( I ) , AVT ( I ) 

C  SVT ( I )  =  FRACTION  OF  PISTON  TRAVEL  (CM). 

C  WILL  BE  NORMALIZED  TO  MAX  PISTON  TRAVEL. 

C  MAX  PISTON  TRAVEL  CHOSEN  TO  AGREE  WITH  INITIAL  LIQUID  VOLUME. 

C  AVT ( I )  =  CORRESPONDING  VENT  AREA  (CM**2). 

C  FIND  VENT  AREA  BY  INTERPOLATING  TABLE. 


81 


oooooooooooooooooooooooonooooooooooooooooooooooooooo 


FILE  8 . NVENT+2 : POPEN 

POPEN  =  LIQUID  PRESSURE  AT  WHICH  VENT  WILL  OPEN. 

IF  VENT 2  THEN  CENTRAL  ROD  OR  BOLT 
WITH  INFINITELY  THIN  PISTON. 

FILE  8.1:  NVENT 

NVENT  =  NO.  OF  TABLE  ENTRIES. 

FILE  8.2. . . :SVT(I) ,AVT(I) 

SVT(I)  =  FRACTION  OF  PISTON  TRAVEL  (CM). 

WILL  BE  NORMALIZED  TO  MAX  PISTON  TRAVEL. 

MAX  PISTON  TRAVEL  CHOSEN  TO  AGREE  WITH  INITIAL  LIQUID  VOLUME. 
AVT(I)  =  CORRESPONDING  VENT  AREA  (CM**2) . 

FILE  8. NVENT+2:  AHOLE , AGRES 

AHOLE  =  AREA  OF  CENTRAL  HOLE  IN  ANNULAR  PISTON  (CM**2) . 

AGRES  =  AREA  OF  GREASE  DYKE. 

INTERPOLATE  CENTRAL  ROD  RADII  TO  FIND  VENT  AREA. 

IF  VENT3  THEN  CENTRAL  ROD  OR  BOLT 
WITH  ACTUAL  PISTON. 

FILE  8.1:  NVENT 

NVENT  =  NO.  OF  TABLE  ENTRIES. 

FILE  8.2...: SVT ( I )  , RROD ( I ) 

SVT(I)  =  DISTANCE  FROM  FRONT  OF  ROD. 

RROD (I)  =  CORRESPONDING  ROD  RADIUS. 

FILE  8. NVENT+2:  AHOLE, AGRES 

AHOLE  =  AREA  OF  CENTRAL  HOLE  IN  ANNULAR  PISTON  (CM* *2). 

AGRES  =  AREA  OF  GREASE  DYKE. 

INITIAL  LIQUID  VOLUME  VI  IS  CONSIDERED  ACCURATE. 

WHEN  THIS  MUCH  LIQUID  HAS  BEEN  INJECTED, 

THE  PISTON  STOPS. 

SOME  PART  OF  THE  PISTON  IS  ASSUMED  TO  HAVE  HIT 
THE  BACK  WALL. 

ACTUAL  SHAPE  OF  PISTON  IS  IRRELEVANT. C 
FILE  8.NVENT+3:  NBLOC 

IF  NBLOC  >  0,  THEN  BLOCK  IS  MOUNTED  ON  BELLEVILLE  SPRINGS. 
BLOCK  WILL  MOVE  BACK  UNTIL  SPRINGS  BOTTOM  OUT. 

FILE  8 . NVENT+4 . . , : SBLOC ( I ) , FBLOC ( I ) 

SBLOC(I)  =  DISTANCE  SPRINGS  ARE  COMPRESSED  (CM). 

FBLOC (I)  =  FORCE  EXERTED  BY  SPRINGS  (MPA-CM**3) . 

FILE  8. NVENT+4 +NBLOC:  ABLOC , BLWT 
ABLOC  =  AREA  OF  BLOCK  (CM**2). 

BLWT  *  MASS  OF  BLOCK  (G) . 

FILE  9:  TPIS  -  PISTON  FRICTION  RESISTANCE  OPTION 

IF  PIS1  THEN  RESISTANCE  IS  A  FUNCTION  OF  PISTON  TRAVEL. 

FILE  9.1:  NPIS 

NPIS  *  NO.  OF  ENTRIES  IN  PISTON  RESISTANCE  TABLE. 

FILE  9.1...:  SFRAC ( I ) , PRS ( I ) 

SFRAC(I)  -  FRACTION  OF  PISTON  TRAVEL. 

WILL  BE  NORMALIZED  TO  MAX  PISTON  TRAVEL. 

PRS (I)  -  PISTON  RESISTANCE  (MPA). 
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PISTON  RESISTANCE  IS  FOUND  BY  INTERPOLATING  TABLE. 

FILE  10:  TDIS  -  DISCHARGE  COEFFICIENT  OPTION  SWITCH 

IF  DIS1  THEN  DC  FUNCTION  OF  PISTON  TRAVEL. 

FILE  10.1:  NDC 

NDC  =  NO.  OF  TABLE  ENTRIES. 

FILE  10.2...:  SFRAC(I) ,DCF(I) 

SFRAC(I)  =  FRACTION  OF  PISTON  TRAVEL. 

WILL  BE  NORMALIZED  TO  MAX  PISTON  TRAVEL. 

DCF (I)  ■  DISCHARGE  COEFFICIENT. 

DC  IS  FOUND  BY  INTERPOLATING  TABLE. 

FILE  11:  TDIS34  -  DISCHARGE  COEFFICIENT  OPTION  SWITCH-TUBE 

IF  DIS1  THEN  DC34  FUNCTION  OF  PROJ  TRAVEL. 

FILE  11.1:  NDC 

NDC  =  NO.  OF  TABLE  ENTRIES. 

FILE  11.2...:  SFRAC(I) , DCF (I) 

SPROJ(I)  =  PROJECTILE  TRAVEL. 

DCF (I)  =  DISCHARGE  COEFFICIENT. 

DC34  IS  FOUND  BY  INTERPOLATING  TABLE. 

FILE  12:  TFLUX  -  MASS  FLUX  OPTION  SWITCH  FOR  ORIFICE 

IF  FLUX1  THEN  STEADY  STATE  FORMULATION. 

FLUID  ENTRANCE  VELOCITY  =  0. 

FILE  12.1:  NVO 

NVO  =  NO.  OF  OPENINGS  IN  PISTON. 

USED  TO  FIND  REYNOLDS  NO. 

IF  FLUX2  THEN  UNSTEADY  MASS  FLUX  FORMULATION. 

TAKE  TIME  DERIVATIVE  OF  MASS  FLUX. 

ASSUME  MASS  FLUX  CONSTANT  WRT  X. 

FILE  12.1:  NVO 

NVO  =  NO.  OF  OPENINGS  IN  PISTON. 

USED  TO  FIND  REYNOLDS  NO. 

IF  FLUX 3  THEN  UNSTEADY  MASS  FLUX  FORMULATION. 

TAKE  TIME  DERIVATIVE  OF  INFLUX  VELOCITY. 

HEISER  FORMULATION. 

ASSUME  MASS  FLUX  CONSTANT  WRT  X. 

FILE  12.1:  NVO, PISTH 

NVO  =  NO.  OF  OPENINGS  IN  PISTON. 

USED  TO  FIND  REYNOLDS  NO. 

PISTH  =  THICKNESS  OF  PISTON  (CM) . 

FILE  13:  TFLUXP  -  MASS  FLUX  OPTION  SWITCH  FOR  TUBE 

IF  FLUX1  THEN  STEADY  STATE  FORMULATION. 

GENERALIZED  BERNOULLI  WITH  ENERGY  DIFFERENCE. 

ISENTROPIC  FLOW. 
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FLUID  ENTRANCE  VELOCITY  =  0. 

IF  FLUX2  THEN  STEADY  STATE  FORMULATION. 

ISENTROPIC  FLOW. 

BUT  SOLVE  FOR  PL  INSTEAD  OF  U34. 

USE  WITH  TUBE3  OR  TUBE4  OPTION. 

FILE  14:  TPROJ  -  PROJECTILE  RESISTANCE  OPTION  SWITCH 

IF  PROJ1  THEN  PROJECTILE  RESISTANCE  FUNCTION  OF  PROJ  TRAVEL. 

FILE  14.1:  NPROJ 

NPROJ  =  NO.  OF  TABLE  ENTRIES. 

FILE  14.2...: STR ( I ) , PTR ( I ) 

STR(I)  =  PROJECTILE  TRAVEL  (CM). 

PTR (I)  =  RESISTIVE  PRESSURE  (MPA). 

INTERPOLATE  TABLE  TO  FIND  PROJECTILE  RESISTANCE  PRESSURE. 
FILE  15:  RH1I , RK1 , RK2 

RH1I  =  LIQUID  DENSITY  AT  ZERO  PRESSURE  (GM/CM**3). 

RK1  =  BULK  MODULUS  AT  ZERO  PRESSURE  (MPA) . 

RK2  =  DERIVATIVE  OF  BULK  MODULUS. 

FILE  16:  ENER , GAM 

ENER  =  CHEMICAL  ENERGY  OF  PROPELLANT  (JOULES/GM) . 

GAM  =  RATIO  OF  SPECIFIC  HEATS. 

FILE  17:  SIGMA, VISK 

SIGMA  =  SURFACE  TENSION  OF  PROPELLANT  (DYNES/ CM) . 

VISK  =  KINEMATIC  VISCOSITY  (CM**2/SEC) . 

FILE  18:  WG , COV 

WG  =  MOLECULAR  WEIGHT  OF  THE  GASES  (GM/MOLE) . 

COV  =  COVOLUME  OF  THE  GAS  (CM**3/GM) . 

FILE  19:  PLI , PGI 

PLI  =  INITIAL  LIQUID  PRESSURE  (MPA) . 

PGI  «  INITIAL  GAS  PRESSURE  (MPA) . 

FILE  20:  TDROP  -  DROPLET  OPTION  SWITCH 

IF  DROP1  THEN  INSTANTANEOUS  BURNING. 

IF  DROP2  THEN  FIXED  SIZE  DROPLETS. 

FILE  20.1: DDR, PBR 

DDR  =  DIAMETER  OF  DROPLETS  (CM) . 

PBR  =  BREAK  PRESSURE  FOR  CHANGE  IN  BURNING  RATE  (MPA) . 

FILE  20.2: ADR1 , BDR1 

BURNING  RATE  =  ADR1  *  (P**DBR1) 

IF  P<PBR . 

FILE  20.3: ADR2 , BDR2 

BURNING  RATE  =  ADR2  *  (P**DBR2) 

IF  P>PBR. 
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c 

C  IF  DROP 3  THEN  FIXED  SIZE  DROPLETS  WRT  SPACE. 

C  DROPLET  SIZE  CHANGES  WITH  SIZE. 

C  FILE  20.1: DDR , PBR 

C  DDR  IS  IGNORED. 

C  PBR  =  BREAK  PRESSURE  FOR  CHANGE  IN  BURNING  RATE  (MPA) . 

C  FILE  20 . 2 : ADR1 , BDR1 

C  BURNING  RATE  =  ADR1  *  (P**DBR1) 

C  IF  P<PBR. 

C  FILE  20.3: ADR2 , BDR2 

C  BURNING  RATE  =  ADR2  *  (P**DBR2) 

C  IF  P>PBR. 

C  FILE  20.4 :NDIAM 

C  NDIAM  =  NO.  OF  TABLE  ENTRIES. 

C  FILE  20.5. .. :SFRAC(I) ,DIAM(I) 

C  SFRAC(I)  =  FRACTION  OF  PISTON  TRAVEL. 

C  WILL  BE  NORMALIZED  TO  MAX  PISTON  TRAVEL. 

C  DIAM(I)  =  DROPLET  DIAMETER 

C  THE  DROPLET  DIAMETER  AT  ANY  TIME  WILL  BE  FOUND  BY 

C  INTERPOLATING  THE  ABOVE  TABLE. 

C 

C  FILE  21:  TPRIM  -  PRIMER  OPTION  SWITCH 

C 

C  IF  PRIM1  THEN  PRIMER  INSTANTANEOUSLY  GOES  TO  GAS. 

C  THE  PRIMER  MASS  IS  THE  MASS  OF  GAS  IN  THE  CHAMBER. 

C 

C  IF  PRIM2  SET  PRIMER  AS  LIQUID  DROPLETS. 

C  ASSUME  PRIMER  HAS  SAME  PROPERTIES  AS  LIQUID  PROPELLANT. 

C  FILE  21.1:  RMPRIM 

C  RMPRIM  =  MASS  OF  PRIMER  (G) . 

C  SET  PRIMER  AS  LIQUID  DROPLETS  IN  COMBUSTION  CHAMBER. 

C  THE  INITIAL  GAS  PRESSURE  IS  ASSUMED  TO  HAVE  BEEN 

C  GENERATED  BY  PRIMER  THAT  HAS  ALREADY  BURNED. 

C 

C  IF  PRIM3  THEN  INJECT  PRIMER  AT  A  CONSTANT  RATE. 

C  ASSUME  PRIMER  HAS  SAME  PROPERTIES  AS  LIQUID  PROPELLANT. 

C  FILE  21.1:  RMPRIM,  TINJEC 

C  RMPRIM  =  MASS  OF  PRIMER  (G) . 

C  TINJEC  =  TIME  FOR  COMPLETE  INJECTION  (S) . 

C 

C  FILE  22:  THEATC  -  HEAT  LOSS  OPTION  (CHAMBER) 

C 

C  IF  HEAT1  THEN  NO  H’TYT  LOSS  TO  THE  CHAMBER  WALLS. 

C 

C  IF  HEAT 2  THEN  CONVECTIVE  HEAT  LOSS  TO  THE  CHAMBER. 

C  Q  «  H  *  A  *  DELTA  T. 

C  FILE  22.1: TEMPC 

C  TEMPC  =  CHAMBER  WALL  TEMPERATURE  (K) . 

C  FILE  22.2 :NCHAM 

C  NCHAM  =  NUMBER  OF  HEAT  LOSS  TERMS. 

C  FILE  22.3. . .25. NCHAM+2 : SPS , HTRAN 

C  SPS  =  PISTON  TRAVEL  (CM) . 
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HTRAN  =  HEAT  TRANSFER  COEFF.  TIMES  AREA  (JOULES/S-K) . 

FILE  2 3 : THEAT  -  HEAT  LOSS  OPTION  (GUN  TUBE) 

IF  HEAT1  THEN  NO  HEAT  LOSS  TO  THE  TUBE  WALLS. 

IF  HEAT 2  THEN  CONVECTIVE  HEAT  LOSS  TO  THE  GUN  TUBE. 

USE  FITS  FOR  VISCOSITY  AND  THERMAL  CONDUCTIVITY. 

FITS  WERE  MADE  FOR  HAN1845  BETWEEN  1500K  AND  3500K. 

FILE  23.1: TEMPW , HLFAC 

TEMPW  =  TUBE  WALL  TEMPERATURE  (K) . 

HLFAC  =  FUDGE  FACTOR. 

FILE  24:  TSHOCK 

IF  SHOCK1  NO  AIR  SHOCK  TERM. 

IF  SHOCK2  INCLUDE  AIR  SHOCK  EFFECT. 

FILE  24.1: AIRP , AIRT 

AIRP  =  INITIAL  PRESSURE  IN  GUN  TUBE  (MPA) . 

AIRT  =  INITIAL  TEMPERATURE  IN  GUN  TUBE  (K) . 

FILE  24.2: AIRGAM , AIRMW 

AIRGAM  =  RATIO  OF  SPECIFIC  HEATS  OF  GAS  IN  BARREL. 

AIRMW  =  MOLECULAR  WEIGHT  OF  GAS  IN  BARREL  (GM/MOLE) . 

FILE  25:  TTUBE 

IF  TUBE1  THEN  STANDARD  LAGRANGE  DISTRIBUTION. 

IF  TUBE2  THEN  LAGRANGE  WITH  NON-ZERO  ENTRANCE  VELOCITY. 

ASSUME  DENSITY  IS  CONSTANT  WRT  X. 

IF  TUBE 3  THEN  UNSTEADY  LAGRANGE  WITH  NON-ZERO  ENTRANCE  VELOCIT 
ASSUME  DENSITY  IS  CONSTANT  WRT  X. 

SET  UP  ODE  FOR  U34. 

USE  WITH  FLUX2  OPTION. 

FILE  25.1:  TTUBE 

IF  BURN1  THEN  TRACK  RAREFACTION  WAVE  AFTER  BURNOUT. 

OTHERWISE  DO  NOT. 

IF  TUBE 4  THEN  UNSTEADY  LAGRANGE  WITH  NON-ZERO  ENTRANCE  VELOCIT 
ASSUME  DENSITY  IS  LINEAR  WRT  X. 

SET  UP  ODE  FOR  U34. 

USE  WITH  FLUX2  OPTION. 

FILE  25.1:  TTUBE 

IF  BURN1  THEN  TRACK  RAREFACTION  WAVE  AFTER  BURNOUT. 

OTHERWISE  DO  NOT. 

FILE  26:  TBUFF 

IF  BUFFI  THEN  NO  WATER  BUFFER. 
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C  IF  BUFF2  THEN  WATER  BUFFER. 

C  FILE  26.1:  APIS5 ,  A0UT5 

C  APIS 5  =  AREA  OF  BACK  OF  PISTON  IN  BUFFER. 

C  AOUT5  =  AREA  OF  EXIT  FROM  BUFFER  TO  OUTSIDE. 

C  FILE  26.2:  CD5 

C  CD5  =  DISCHARGE  COEFFICIENT  INTO  OUTSIDE. 

C  FILE  26.3:  V5 

C  V5  =  INITIAL  VOLUME  OF  THE  BUFFER. 

C  FILE  26.4:  P5 

C  P5  =  INITIAL  PRESSURE  OF  THE  BUFFER. 

C  FILE  26.5:  RH50 ,  RK51,  RK52 

C  RH50  =  DENSITY  AT  ZERO  PRESSURE. 

C  RK51  *  BULK  MODULUS  AT  ZERO  PRESSURE. 

C  RK52  =  DERIVATIVE  OF  BULK  MODULUS. 

C 

C  FILE  27:  TINC,HTOP 

C  TINC  =  TIME  INCREMENT  BETWEEN  OUTPUT  LINES  (SEC) . 

C  HTOP  =  MAX  TIME  STEP  ALLOWED  IN  THE  INTEGRATION  (SEC) . 

C 

C  FILE  28:  EPS,SREC 

C  EPS  =  ERROR  CONTROL  FOR  TIME  INTEGRATION. 

C  SREC  =  CUTOFF  BETWEEN  RELATIVE  AND  ABSOLUTE  ERROR  CONTROL. 

C 

C  FILE  2 9 : MF , KWRITE 

C  MF  -  METHOD  OF  INTEGRATION. 

C  USUALLY  MF=22  (STIFF  INTEGRATION  SCHEME  WITH  INTERNALLY 

C  GENERATED  JACOBIAN) . 

C  KWRITE  ■  DIAGNOSTIC 

C  IF  KWRITE=1,  THEN  PRINT  OUT  AFTER  EACH  TIME  STEP. 

C  ALSO  TURN  ON  EPISODE  ERROR  MESSAGES. 

C 

C  FILE  30 : TMAX 

C  TMAX  =  MAX  RUN  TIME  ALLOWED. 

C 

C  FILE  31 : TREP 

C 

C  IF  REP1  THEN  INTEGRATE  ONCE. 

C 

C  IF  REP2  ITERATE:  CHANGE  VENT  AREA  TO  ADJUST  LIQUID  PRESSURE. 

C  FILE  31.1:  PTAR,  VINC 

C  PTAR  =  DESIRED  LIQUID  PRESSURE  (MPA) . 

C  VINC  =  INITIAL  INCREMENT  FOR  CHANGING  VENT  AREA  (CM**2). 

C 

C  IF  REP3  ITERATE:  CHANGE  RESERVOIR  VOLUME  TO  ADJUST  MUZZLE  V. 

C  FILE  31.1:  VMTAR,  VINC 

C  VMTAR  =  DESIRED  MUZZLE  VELOCITY  (M/S) . 

C  VINC  =  INITIAL  INCREMENT  FOR  CHANGING  RESERVOIR  VOLUME  (CM* *3 

C 

C  FILE  3 2 : TCHAM 

C 

C  IF  CHAM1  THEN  COMPUTE  CHAMBER  PRESSURE. 

C 
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IF  CHAM2  THEN  USE  EXPERIMENTAL  CHAMBER  PRESSURE. 
READ  IN  FROM  STANDARD  GRAPHICS  FILE  (TAPE1) . 
FIRST  COLUMN  =  TIME  (MS). 

THIRD  COLUMN  -  CHAMBER  PRESSURE  (MPA) . 
*********************************************** 
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APPENDIX  B 


Below  is  a  listing  of  the  job  stream  for  the  first  test  problem 
(instantaneous  burning  with  constant  discharge  coefficient) . 
Following  is  a  listing  of  the  output.  Since  there  are  many 
variables  to  keep  track  of,  the  output  is  on  multiple  files. 
Only  the  ones  of  interest  for  a  particular  problem  need  to  be 
printed  out.  Here  only  the  minimal  output  is  given. 


30LC21  -  2/3 
0.0 
3.0 
287.1 
2109.1 
166.98 
33.778 


Charge  - 
243.84 


95.0 

45.508 


Instan  -  CD=0.95 

OFFSET  PROJ  TRAVEL 
GUN  TUBE  DIAMETER 
PROJ  WEIGHT 
PISTON  WEIGHT 
VI  V3 
A1  A3 


VENT  3 

CENTRAL  BOLT  WITH 

7 

BOLT  RADIUS  VER! 

0.0 

1.826 

0.127 

1.826 

0.363 

1.800 

0.558 

1.800 

1.320 

1.651 

5.889 

1.651 

8.103 

1.822 

10.475 

.666 

AHOLE  AGRES 

11 

NBLOC 

.000 

163.443 

SBLOC  FBLOC 

.007 

171.470 

.034 

202.162 

.058 

229.517 

.085 

265.323 

.  110 

288.898 

.135 

316.253 

.  160 

344.275 

.  185 

370.963 

.210 

398.318 

.238 

427.675 

23.349 

3921.0 

ABLOC  BLOCK  1 

PIS1 

2 

PISTON  RESISTANCE 

0.0 

0.0 

1.0 

0.0 

DIS1 

2 


0.0 

1.0 


.95 

.95 


DIS1 

2 


DIS.  COEFF.  VERSUS  PISTON  TRAVEL 


DIS.  COEFF.  VS.  PROJ  TRAVEL  -  TUBE 


0.0 


1.0 
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243.84  1.0 

FLUX1 

1  1.0 
FLUX1 
PROJ1 
3 

0.0  50.0 

0.01  5.0 

243.84  5.0 

1.43  5350.0  9.11 


4035.5 
66.9 

22.848 

7.0 

DROP1 

PRIM3 

1.5  0 
HEAT1 
HEAT2 

300.0 
SHOCK2 
0.1 
1.4 
TUBE  2 
BUFFI 

1.00E-04 
1.00E-06 
22  0 
40.0 

REP1 

CHAM1 


1.2226 

.04988 

.677 

0.1 


0025 


1.0 

300. 

28.84 


1.00E-05 

1.00E-09 


STEADY  STATE  MASS  FLUX  FORMULATION 
NVO  PTH 

ISENTROPIC  FLOW  INTO  TUBE 
PROJ  RESISTANCE 


RHO  K1  K2 
ENERGY  GAMMA 

SURFACE  TENSION  KINEMATIC  VISCOSITY 
MOL  WT  GAS  COVOLUME 

PI  P3 

INSTANTANEOUS  BURNING 
INJECT  HOT  GAS 

PRIMER  MASS  INJECTION  TIME 
NO  HEAT  LOSS  TO  CHAMBER  WALLS 
HEAT  LOSS  TO  GUN  TUBE  WALLS 
TUBE  TEMP  FUDGE  FACTOR 
AIR  SHOCK 
AIRP  AIRT 
AIRGAM  AIRMW 

MODIFIED  LAGRANGE  DISTRIBUTION 

NO  WATER  BUFFER 

TINC  HTOP 

EPS  SREC 

MF  KWRITE 

TMAX 

INTEGRATE  ONCE 
COMPUTE  P3 
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P  RO  J  1  PRO J  RESISTANCE 


EPS  *  1 . OOOOOE -  06  SPEC  =  1.00000E-09 
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MUZZLE  VEL  (M/SEC)  1083 
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APPENDIX  C 


Below  is 

a  listing  of 

the  job  stream  for  the  second  test 

problem 

(droplet 

formation  with  constant  discharge  coefficient) . 

Following 

is  the  summary  page  from  the  output. 

30LC31  - 

2/3  Charge  - 

Drop  -  CD=0.95 

0.0 

243.84 

OFFSET  PROJ  TRAVEL 

3.0 

GUN  TUBE  DIAMETER 

287.1 

PROJ  WEIGHT 

2109.1 

PISTON  WEIGHT 

166.98 

95.0 

VI  V3 

33.778 

45.508 

A1  A3 

VENT  3 

CENTRAL  BOLT  WITH  ACTUAL  PISTON 

7 

BOLT  RADIUS  VERSUS  PISTON 

TRAVEL 

0.0 

1.826 

0.127 

1.826 

0.363 

1.800 

0.558 

1.800 

1.320 

1.651 

5.889 

1.651 

8.103 

1.822 

10.475 

.  666 

AHOLE  AGRES 

11 

NBLOC 

.  000 

163.443 

SBLOC  F3L0C 

.007 

171.470 

.034 

202 . 162 

.058 

229.517 

.085 

265.323 

.  110 

288.898 

.135 

316.253 

.  160 

344.275 

.  185 

370.963 

.210 

398.318 

.238 

427 . 675 

23.349  3921.0 

ABLOC  BLOCK  WEIGHT 

PIS1 

O 

PISTON  RESISTANCE 

0.0 

0.0 

1.0 

0.0 

DIS1 

o 

DIS.  COEFF.  VERSUS  PISTON  TRAVEL 

0.0 

.95 

1.0 

.95 

DIS1 

DIS.  COEFF.  VS.  PROJ  TRAVEL 

-  TUBE 

£> 

0.0 

1.0 

243.84 

1.0 

FLUX1 

STEADY  STATE  MASS  FLUX  FORMULATION 

1 

1.0 

NVO  PTH 

99 


FLUX1 

ISENTROPIC  FLOW  INTO  TUBE 

PROJ1 

"i 

PROJ  RESISTANCE 

J 

0.0 

50.0 

0.01 

5.0 

243.84 

5.0 

1.43  5350.0  9.11 

RHO  K1  K2 

4035.5 

1.2226 

ENERGY  GAMMA 

66.9 

.04988 

SURFACE  TENSION  KINEMATIC  VISCOSITY 

22.848 

.677 

MOL  WT  GAS  COVOLUME 

7.0 

0.1 

PI  P3 

DROP  3 

MONOSIZE  DROPS-  FUNC  OF  PISTON  TRAVEL 

O.Ol 

95.2590 

DDR  PBR 

1.64 

.  103 

ADR1  BDR1 

1.64 

*7 

.  103 

ADR1  BDR1 

/ 

0.00 

1.00 

0.363 

1.00 

0.364^ 

0.04 

0.558  — 

0.02 

1.320 

0.016 

3.000 

0.008 

6.9227 

0.008 

PRIM3 

INJECT  HOT  GAS 

1.5  0. 

0025 

PRIMER  MASS  INJECTION  TIME 

HEAT1 

NO  HEAT  LOSS  TO  CHAMBER  WALLS 

HEAT2 

HEAT  LOSS  TO  GUN  TUBE  WALLS 

300.0 

1.0 

TUBE  TEMP  FUDGE  FACTOR 

SHOCK2 

AIR  SHOCK 

0.1 

300. 

AIRP  AIRT 

1.4 

28.84 

AIRGAM  AIRMW 

TUBE2 

MODIFIED  LAGRANGE  DISTRIBUTION 

BUFFI 

NO  WATER  BUFFER 

1.00E-04 

1.00E-05 

TINC  HTOP 

1.00E-06 

1.00E-09 

EPS  SREC 

22  0 

MF  KWRITE 

40.0 

TMAX 

REP1 

INTEGRATE  ONCE 

CHAM1 

COMPUTE  P3 
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MUZZLE  VEL  (M/SEC)  1050.4 

MAX  V  PIS  (M/SEC)  24.1 

MAX  PI  (MPA)  266.4 

MAX  P3  (MPA)  176.2 

MAX  PL  (MPA)  159.1 

MAX  PR  (MPA)  150.6 

MAX  ACC  (K-G)  33.9 

MAX  MASS  ERROR  0.00 

MAX  ENERGY  ERROR  0.18 

BALLISTIC  EFFICIENCY  =  16.42  % 

EXPANSION  RATIO  =  7.91 


LOSS  TO  CHAMBER  WALLS  =  11.49  % 


RUN  TIME  =  5.6  NSTEP  = 


1496 
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APPENDIX  D 


Below  is  a  listing  of  the  job  stream  for  the  third  test  problem 
(instantaneous  burning;  constant  discharge  coefficient;  longer 
Belleville) .  Following  is  the  summary  page  from  the  output. 


30LC22  - 

2/3  Charge  - 

Instan  -  CB=0.95  -  Long  Belle. ille 

0.0 

243.84 

OFFSET  PROJ  TRAVEL 

3.0 

GUN  TUBE  DIAMETER 

287.1 

PROJ  WEIGHT 

2109.1 

PISTON  WEIGHT 

166.98 

95.0 

VI  V3 

33.778 

45.508 

A1  A3 

VENT  3 

CENTRAL  BOLT  WITH  ACTUAL  PISTON 

7 

BOLT  RADIUS  VERSUS  PISTON  TRAVEL 

0.0 

1.826 

0.127 

1.826 

0.363 

1.800 

0.558 

1.800 

1.320 

1.651 

5.889 

1.651 

8.103 

1.822 

10.475 

.666 

AHOLE  AGRES 

12 

NBLOC 

.000 

163.443 

SBLOC  FBLOC 

.007 

171.470 

.034 

202.162 

.058 

229,517 

.085 

265.323 

.  110 

288.898 

.135 

316.253 

.160 

344.275 

.  185 

370.963 

.210 

398.318 

.238 

427.675 

.438 

500.000 

23,349  3921.0 

ABLOC  BLOCK  WEIGHT 

PIS1 

n 

PISTON  RESISTANCE 

0.0 

0.0 

1.0 

0.0 

DIS1 

DIS.  COEFF.  VERSUS  PISTON  TRAVEL 

& 

0.0 

.95 

1.0 

.95 

DIS1 

o 

DIS.  COEFF.  VS.  PROJ  TRAVEL  -  TUBE 

€» 

0.0 

1.0 

243.84 

1.0 

FLUX1 

STEADY  STATE  MASS  FLUX  FORMULATION 
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1.0 


1 

FLUX1 

PROJ1 

3 


0.0  50.0 


0.01  5.0 

243.84  5.0 

1.43  5350.0  9.11 


4035.5 
66.9 

22.848 

7.0 

DROP1 

PRIM3 

1.5  0 
HEAT1 
HEAT  2 

300.0 
SHOCK2 
0.1 
1.4 
TUBE  2 
BUFFI 

1 . 00E-04 
1.00E-06 
22  0 
40.0 

REP1 

CHAM1 


1.2226 

.04988 

.677 

0.1 


0025 


1.0 

300. 

28.84 


1.00E-05 

1.00E-09 


NVO  PTH 

ISENTROPIC  FLOW  INTO  TUBE 
PROJ  RESISTANCE 


RHO  K1  K2 
ENERGY  GAMMA 

SURFACE  TENSION  KINEMATIC  VISCOSITY 
MOL  WT  GAS  COVOLUME 

PI  P3 

INSTANTANEOUS  BURNING 
INJECT  HOT  GAS 

PRIMER  MASS  INJECTION  TIME 
NO  HEAT  LOSS  TO  CHAMBER  WALLS 
HEAT  LOSS  TO  GUN  TUBE  WALLS 
TUBE  TEMP  FUDGE  FACTOR 
AIR  SHOCK 
AIRP  AIRT 
AIRGAM  AIRMW 

MODIFIED  LAGRANGE  DISTRIBUTION 

NO  WATER  BUFFER 

TINC  HTOP 

EPS  SREC 

MF  KWRITE 

TMAX 

INTEGRATE  ONCE 
COMPUTE  P3 
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MUZZLE  VEL  (M/SEC)  1140.9 

MAX  V  PIS  (M/SEC)  28.0 

MAX  PI  (MPA)  351.8 

MAX  P3  (MPA)  238.1 

MAX  PL  (MPA)  214.5 

MAX  PR  (MPA)  210.8 

MAX  ACC  (K-G)  50.0 

MAX  MASS  ERROR  0.00 

MAX  ENERGY  ERROR  0.14 


BALLISTIC  EFFICIENCY  =  19.37  % 

EXPANSION  RATIO  =  7.94 

LOSS  TO  CHAMBER  WALLS  =  11.08  % 


RUN  TIME  =  3.8  NSTEP  = 


1169 
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APPENDIX  E 


Below  is 

a  listing  of  the 

job  stream  for  the  fourth  test  problem 

(droplet 

formation;  constant  discharge  coefficient;  longer 

Belleville) .  Following  is 

the  summary  page  from  the  output. 

3  0  LC  32- 

2/3  Charge  -  Drop 

-  CD=0 . 95  -  Longer  Belleville 

0.0 

243.84 

OFFSET  PROJ  TRAVEL 

3.0 

GUN  TUBE  DIAMETER 

287.1 

PROJ  WEIGHT 

2109.1 

PISTON  WEIGHT 

166.98  95.0 

VI  V3 

33.778  45.508 

A1  A3 

VENT  3 

CENTRAL  BOLT  WITH  ACTUAL  PISTON 

7 

BOLT  RADIUS  VERSUS  PISTON  TRAVEL 

0.0 

1.826 

0.127 

1.826 

0.363 

1.800 

0.558 

1.800 

1.320 

1.651 

5.889 

1.651 

8.103 

1.822 

10.475 

.666 

AHOLE  AGRES 

12 

NBLOC 

.000 

163.443 

SBLOC  FBLOC 

.007 

171.470 

.034 

202.162 

.058 

229.517 

.085 

265.323 

.110 

288.898 

.135 

316.253 

.160 

344.275 

.185 

370.963 

.210 

398.318 

.238 

427.675 

.438 

500. 

23.349 

3921.0 

ABLOC  BLOCK  WEIGHT 

PIS1 

O 

PISTON  RESISTANCE 

Z 

0.0 

0.0 

1.0 

0.0 

DIS1 

DIS.  COEFF.  VERSUS  PISTON  TRAVEL 

£ 

0.0 

.95 

1.0 

.95 

DIS1 

DIS.  COEFF.  VS.  PROJ  TRAVEL  -  TUBE 

z 

0.0 

1.0 

243.84 

1.0 

FLUX1 

STEADY  STATE  MASS  FLUX  FORMULATION 

1.0 


1 

FLUX1 

PROJ1 

3 

0.0 


50.0 


0.01 

5.0 

243.84 

5.0 

1.43  5350.0  9.11 

4035.5 

1.2226 

66.9 

.04988 

22.848 

.677 

7.0 

0.1 

DROP3 

0.01 

95.2590 

1.64 

.103 

1.64 

.  103 

7 

0.00 

1.00 

0.558 

1.00 

0.559 

0.05 

1.320 

0.05 

3.000 

0.02 

5.889 

0.012 

7.115 

0.01 

PRIM3 

1.5  0. 

0025 

HEAT1 

HEAT  2 

300.0 

1.0 

SHOCK2 

0.1 

300. 

1.4 

28.84 

TUBE2 

BUFFI 

1.00E-04 

1.00E-05 

1.00E-06 

1.00E-09 

22  0 

40.0 

REP1 

CHAM1 

NVO  PTH 

ISENTROPIC  FLOW  INTO  TUBE 
PROJ  RESISTANCE 


RHO  K1  K2 
ENERGY  GAMMA 

SURFACE  TENSION  KINEMATIC  VISCOSITY 
MOL  WT  GAS  COVOLUME 

PI  P3 

MONOSIZE  DROPS-  FUNC  OF  PISTON  TRAVEL 
DDR  PBR 

ADR1  BDR1 

ADR1  BDR1 


INJECT  HOT  GAS 

PRIMER  MASS  INJECTION  TIME 
NO  HEAT  LOSS  TO  CHAMBER  WALLS 
HEAT  LOSS  TO  GUN  TUBE  WALLS 
TUBE  TEMP  FUDGE  FACTOR 
AIR  SHOCK 
AIRP  AIRT 
AIRGAM  AIRMW 

MODIFIED  LAGRANGE  DISTRIBUTION 

NO  WATER  BUFFER 

TINC  HTOP 

EPS  SREC 

MF  KWRITE 

TMAX 

INTEGRATE  ONCE 
COMPUTE  P3 
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MUZZLE  VEL  (M/SEC)  1058.8 

MAX  V  PIS  (M/SEC)  24.2 

MAX  PI  (MPA)  288.1 

MAX  P3  (MPA)  181.6 

MAX  PL  (MPA)  160.1 

MAX  PR  (MPA)  147.6 

MAX  ACC  (K-G)  32.9 

MAX  MASS  ERROR  0.00 

MAX  ENERGY  ERROR  0.22 

BALLISTIC  EFFICIENCY  =  16.68  % 

EXPANSION  RATIO  =  7.94 


LOSS  TO  CHAMBER  WALLS  =  11.63  % 


RUN  TIME  =  5.4  NSTEP  = 


1485 
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APPENDIX  F 


0.558 

0.800 

1.320 

7.115 


,95 
,  30 
95 
,95 


Below  is 

a  listing  of  the 

job  stream  for  the  fifth  test  problem 

(droplet 

formation;  adjusted  discharge  coefficient;  longer 

Belleville) .  Following  is 

the  summary  page  from  the  output. 

30LC42- 

2/3  Charge  -  Drop 

-  Simple  CD  -  Long  Belle 

0.0 

243.84 

OFFSET  PROJ  TRAVEL 

3.0 

GUN  TUBE  DIAMETER 

287.1 

PROJ  WEIGHT 

2109.1 

PISTON  WEIGHT 

166.98  95.0 

VI  V3 

33.778  45.508 

A1  A3 

VENT3 

CENTRAL  BOLT  WITH  ACTUAL  PISTON 

7 

BOLT  RADIUS  VERSUS  PISTON  TRAVEL 

0.0 

1.826 

0.127 

1.826 

0.363 

1.800 

0.558 

1.800 

1.320 

1.651 

5.889 

1.651 

8.103 

1.822 

10.475 

.666 

AHOLE  AGRES 

12 

NBLOC 

.000 

163.443 

SBLOC  FBLOC 

.007 

171.470 

.034 

202.162 

.058 

229.517 

.085 

265.323 

.110 

288.898 

.135 

316.253 

.160 

344.275 

.185 

370.963 

.210 

398.318 

.238 

427.675 

.438 

500. 

23.349 

3921.0 

ABLOC  BLOCK  WEIGHT 

PIS1 

O 

PISTON  RESISTANCE 

0.0 

0.0 

1.0 

0.0 

DIS1 

c 

DIS.  COEFF.  VERSUS  PISTON  TRAVEL 

0.0 

.95 

DIS1 

2 


DIS.  COEFF.  VS.  PROJ  aKaVEL  -  TUBE 
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0.0 

1.0 

243.84 

FLUX1 

1.0 

STEADY  STATE  MASS  FLUX  FORMULATION 

1 

1.0 

NVO  PTH 

FLUX1 

ISENTROPIC  FLOW  INTO  TUBE 

PROJ1 

PROJ  RESISTANCE 

J 

0.0 

50.0 

0.01 

5.0 

243.84 

5.0 

1.43 

5350.0  9.11 

RHO  K1  K2 

4035.5 

1.2226 

ENERGY  GAMMA 

66.9 

.04988 

SURFACE  TENSION  KINEMATIC  VISCOSITY 

22.848 

.677 

MOL  WT  GAS  CC VOLUME 

7.0 

0.1 

PI  P3 

DROP3 

MONOSIZE  DROPS-  FUNC  OF  PISTON  TRAVEL 

0.01 

95.2590 

DDR  PBR 

1.64 

.103 

ADR1  BDR1 

1.64 

8 

0.00 

.  103 

ADR1  BDR1 

1.0000 

0.558 

1.0000 

0.559 

0.0500 

0.800 

0.0300 

1.320 

0.0150 

3.000 

0.0100 

5.889 

0.0075 

7.115 

PRIM3 

0.0075 

INJECT  HOT  GAS 

1.5 

0.0025 

PRIMER  MASS  INJECTION  TIME 

HEAT1 

NO  HEAT  LOSS  TO  CHAMBER  WALLS 

HEAT  2 

HEAT  LOSS  TO  GUN  TUBE  WALLS 

300.0 

1.0 

TUBE  TEMP  FUDGE  FACTOR 

SHOCK2 

AIR  SHOCK 

0.1 

300. 

AIRP  AIRT 

1.4 

28.84 

AIRGAM  AIRMW 

TUBE 2  MODIFIED  LAGRANGE  DISTRIBUTION 

BUFFI  NO  WATER  BUFFER 

1 . 00E-04  1.00E-05  TINC  HTOP 

1.00E-06  1.00E-09  EPS  SREC 

22  0  MF  KWRITE 

40.0  TMAX 

REP1  INTEGRATE  ONCE 

CHAM1  COMPUTE  P3 
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MUZZLE  VEL  (M/SEC)  1059.5 

MAX  V  PIS  (M/SEC)  24.2 

MAX  PI  (MPA)  277.7 

MAX  P3  (MPA)  177.9 

MAX  PL  (MPA)  158.5 

MAX  PH  (MPA)  149.8 

MAX  ACC  (K-G)  33.8 

MAX  MASS  ERPOR  0.00 

MAX  ENERGY  ERROR  0.15 

BALLISTIC  EFFICIENCY  =  16.70  % 

EXPANSION  RATIO  «  7.94 


LOSS  TO  CHAMBER  WALLS  =  11.47  % 


RUN  TIME  =  5.4  NSTEP  = 


1 


1471 
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This  laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  It  publishes.  Your 
comments /answers  below  will  aid  us  in  our  efforts. 

1.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for  which 
the  report  will  be  used.)  _ 


2.  How,  specifically,  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas, 
etc.) _ _ _ 


3.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaoorate.  _ _ _ 


4.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to 
organ! zation,  technical  content,  format,  etc.)  _ _ _ _ 
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